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1 .  INTRODUCTION 

Let  > ("is  preferred  to")  be  a binary  relation  on  the  set  of  simple 
probability  measures  or  'gambles'  defined  on  a set  T of  consequences. 
Throughout  this  study  it  will  be  assumed  that: 

1.  T is  the  Cartesian  product  of  two  or  more  nondegenevate  closed  and 
bounded  real  intervals; 

2.  > on  & satisfies  the  axioms  of  von  Neumann  and  Morgenstern  (1947) 

or  an  equivalent  system  (Herstein  and  Milnor,  195);  Fishburn,  1970)  so  that 
there  exists  u:  T Re  such  that 

P>Q  iff  £P(t)u(t)  > EQ(t)u(t),  for  all  P,Q  € $>  , (1) 

T T 

with  u*  satisfying  (1)  when  u does  iff  u*  is  a positive  affine  transformation 
of  u of  the  form  u*  * ua^  where  u3^(t)  = au(t)  + b,  a > 0; 

3.  The  von  Neumann-Morgenstern  utility  function  u in  (1)  is  continuous 
(Grandmont,  1972;  Foldes,  1972)  in  the  relative  usual  product  topology  for  T. 

The  purpose  of  the  study  is  to  analyze  methods  for  approximating  u.  The 

present  paper  deals  with  two-dimensional  consequence  spaces  T = X * Y;  a 

sequel  paper  will  exc-.mli  •:>.  T = T * T T with  n > 3.  Although  con- 

12  n — 

tinuity  of  u is  not  implied  by  the  axioms  that  are  necessaiy  and  sufficient 
for  (I),  and  plausible  examples  of  discontinuities  are  easy  to  imagine, 
continuity  is  crucial  to  mtv:*-  work  in  approximation  theory  (e.g.  Cheney, 

1966,  and  Lo-entz,  1966)  and  will  be  assumed  here. 
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It  is  one  thing  to  postulate  the  existence  of  continuous  u on  T that 

sfo 

satisfies  (1)  and  is  unique  up  to  transformations  u with  a > 0,  but  quite 
another  thing  to  estimate  u in  an  actual  decision  situation.  Consequently, 
much  effort-  has  been  devoted  to  theoretical  and  methodological  aspects  of 
utility  function  assessment.  The  assessment  and  approximation  of  single- 
variable utility  functions  is  discussed  by  Meyer  and  Pratt  (1968).,  Bradley 
and  Frey  (1975),  and  Ohlson  and  Kallio  (1975)  among  others.  Theoretical 
work  in  multiattribute  or  multivariate  utility  (Fishburn,  1965,  1974; 

Poliak,  1967;  Raiffa,  1969;  Keeney,  1971,  1972a;  Fishburn  and  Keeney,  1974; 

B'arquhar,  1975)  has  focused  on  axioms  for  > on  that  allow  u(t  ,...,t  ) 

i n 

to  be  written  as  a combination  of  functions  defined  on  fewer  than  n attributes, 

such  as  u(t  ,...,t  ) ■ u (t  ) + u (t  ) +...+  u (t  ) or  u(t  ,t  ,...,t  ) = 

12  n ii  22  nn  i 2 n 

u (t  )u  (t  ) . . .u  (t  ) . A desire  to  simplify  the  task  of  utility  assessment 
l l 2 2 n n 

has  motivated  much  of  this  work.  Examples  of  its  application  to  specific 
situations  are  given  by  Raiffa  (1969),  Keeney  (1972b,  1973)  and  Keeney  and 
Nair  (1974).  Sicherman  (1175)  has  developed  an  interactive  computer  program 
for  assessment  in  the  additive  and  multiplicative  cases. 

Research  workers  who  have  been  involved  in  the  development  of  special 
forms  for  multiattribute  utility  functions  realize  of  course  that  the 
independence  axioms  that  characterize  the  special  forms  may  fail  to  hold  in 
a given  situation  because  of  evaluative  interdependencies  among  the  attri- 
butes. Consequently,  there  is  a need  to  explore  the  general  problem  of 
multiattribute  utility  assessment  in  the  absence  of  simplifying  independence 
assumptions.  Although  there  are  several  approaches  to  this  problem,  the 
present  study  will  focus  on  approximations  of  u that  are  written  as  finite 
sums  of  products  of  functions  on  the  individual  attributes.  In  the  context 
of  T = X > Y,  the  form  of  approximation  that  will  be  used  here  is 
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m 

v(x,y)  »■  £ f.(x;g, (y).  (2) 

i=i  1 1 

The  only  direct  assessment  of  u that  will  be  required  for  (2)  involves  either 
the  evaluation  of  u at  a finite  number  of  points  in  T or  th^  evaluation  of  a 
finite  number  of  single-variable  conditional  utility  functions  of  the  form 
u(x»Yj)  and  u(x^,y),  where  x^  is  a fixed  element  in  X and  y^  is  a fixed 
element  in  Y.  The  functions  f^  and  g^  may  involve  the  conditional  utility 
functions  or  they  may  be  specified  independently  of  any  utility  assessment. 

In  later  sections  it  will  be  assumed — as  a first  approximation — that  there 
is  no  error  in  the  assessment  of  the  u values  used  in  the  right  hand  side 
of  (2). 

There  are  three  main  reasons  for  using  (2).  First,  it  is  generally 
conceded  that  it  is  much  easier  to  assess  single-attribute  utility  funccions 
than  to  assess  two-attribute  functions  in  their  full  generality.  Secondly, 
the  sum-ef-products  form  of  (2)  is  computationally  attractive  in  the  context 
of  optimization  algorithms.  Finally,  the  right  hand  side  of  (2)  subsumes 
the  special  forms  of  u(x,y)  that  have  thus  far  been  characterized  by  indepen- 
dence axioms  (Fishburn,  1974).  The  most  general  of  these  is  u(x,y)  = 
u (x)  + u (y)  + h (x)h  (y),  which  agrees  with  (2)  when  m = 3 and  f ■ u , 

12  12  1 j 

f : 1,  £ =h,g  : 1,  g =u  and  g = h . 

2 3 1 1 2 2 3 2 

The  present  paper  is  organized  as  follows.  Some  basic  ideas  from 
approximation  theory  are  briefly  noted  in  the  next  section,  and  a result 
from  this  theory  is  presented  in  our  utility  context.  Tne  simple  additive 
and  multiplicative  approximations  are  examined  in  section  3.  Section  4 then 
discusses  some  elementary  interpolation  approximations  based  on  finite  sets 
of  single-variable  conditional  utility  functions.  Tl,  final  section  examines 
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approximations  that  are  exact  on  a grid  in  X x Y:  that  is,  they  give 

v(xiSy)  = u^.y)  for  every  xi  in  a finite  subset  of  X and  all  y 6 Y,  along 
with  v(x,yj)  = u(x,y_j)  for  every  y^  in  a finite  .subset  of  Y and  all  x 6 X. 

Readers  who  are  familiar  with  the  diversity  and  scope  of  approximation 
theory  will  realize  that  the  present  study  represents  a very  modest  step  in 
the  development  of  a theory  and  methodology  for  the  approximation  and 
assessment  of  multiattribute  utility  functions.  It  is  hoped  that  the  study 
will  elicit  additional  interest  in  the  topic. 

2.  REMARKS  ON  APPROXIMATION  THEORY 

This  section  outlines  a few  basic  ideas  of  approximation  theory, 
comments  on  aspects  of  (2)  that  will  play  a role  throughout  the  paper,  and 
provides  an  application  of  approximation  theory  to  our  utility  context.  A 
broader  introduction  to  approximation  theory  can  be  obtained  from  the  paper 
by  Buck  (1959)  and  the  books  by  Lorentz  (1966)  and  Cheney  (1966).  Other 
suggested  wortts  include  the  papers  by  Rivlin  and  Shapiro  (1961),  Lorentz 
(1972)  and  Jerome  (1973),  the  collections  edited  by  Langer  (1959)  and 
Lorentz  (1973),  and  various  articles  in  the  Journal  of  Approximation  Theory. 

Appropriate  to  our  purposes  let  S = [ 0 , 1 ] n be  the  n-dimensional  unit 
cube  and  lot  C(S)  be  the  real  linear  space  (Kelley  and  Namioka,  1963) 
of  all  continuous  real  valued  functions  on  S.  The  most  commonly  used  norm 
in  approximation  theory  for  measuring  distances  between  functions  in  C(S) 
is  the  uniform  norm 

| |f | I * sup  | f (s ) [ = max  |f(s)|, 
s£S  S 


f € CJ(S), 


(3) 
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Approximations  that  are  exact  on  a grid  in  X x Y:  that  is,  they  give 

- u(x^,y)  for  every  x^  in  a finite  subset  of  X and  all  y € Y,  along 
with  v(x,yj)  ■ u(x,y^)  for  every  y ^ in  a finite  subset  of  Y and  all  x € X. 

Readers  who  are  familiar  with  the  diversity  and  scope  of  approximation 
theory  will  realise  that  the  present  study  represents  a very  modest  step  in 
the  development  of  a theory  and  methodology  for  the  approximation  and 
assessment  of  multiattribute  utility  functions.  It  is  hoped  that  the  study 
will  elicit  additional  interest  in  the  topic. 

2.  REMARKS  ON  APPROXIMATION  THEORY 

This  section  outlines  a few  basic  ideas  of  approximation  theory, 
coaoents  on  aspects  of  (2)  that  will  play  a role  throughout  the  paper,  and 
provides  an  application  of  approximation  theory  to  our  utility  context.  A 
broader  introduction  to  approximation  theory  can  be  obtained  from  the  paper 
by  Buck  (1959)  and  the  books  by  Lorents  (1566)  and  Cheney  (1966).  Other 
suggested  works  include  the  papers  by  Rlvlin  and  Shapiro  (1961),  Lorents 
(1972)  and  Jerome  (1973),  the  collections  edited  by  Langer  (1959)  and 

V 

Lorents  (1973),  and  various  articles  in  the  Journal  of  Approximation  Theory. 

Appropriate  to  our  purposes  let  S ■ [0,l]n  be  the  n-dimenslonal  unit 
cube  and  lot  C(S)  be  the  real  linear  space  (Kelley  and  Namioka,  1963) 
of  all  continuous  real  valued  functions  on  S.  The  most  commonly  used  norm 
in  approximation  theory  for  measuring  distances  between  functions  in  C(S) 
is  the  uniform  nora 


1 1 f 1 1 * sup  | f (s) | - max  | f (s) | , f € C(S), 

3gS  S 


(3) 


** 
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the  latter  equality  arising  from  compactness  and  continuity.  (Least  squares 

minimization  problems  use  the  quadratic  or  Euclidean  norm.)  The  uniform  norm 

will  be  used  throughout  our  study.  Convergence  in  this  norm,  i.e.  j | f m - f| | 

+ 0 as  m + “ with  f,f  ,f  , ...  6 C(S),  is  equivalent  to  the  uniform  convergence 

of  the  f to  f. 
m 

Let  f € C(S)  be  given  and  let  D be  a nonempty  subset  of  C(S)  whose 

functions  are  proposed  as  approximations  of  f.  Define  d(D,f)  * inf  ||g  - f|| 

g€D 

and  let  * {g  € D:  ||g  - f||  » d(D,f)}.  Hence  d(D,f)  is  the  ’distance' 

from  D to  f and  the  functions  in  provide  the  best  approximations  to  f as 
g^’jg^d  by  (3).  Among  otner  things,  approximation  theory  is  concerned  with 
the  following  questions: 

1.  Can  d(D,f)  be  specified  precisely,  or  is  it  possible  to  obtain 
tight  bounds  on  it? 

2.  What  can  be  said  about  the  functions  in  D for  which  | |g  - f | j < 
d(D,f)  + 6 for  given  <5  > 0? 

3.  Is  nonempty,  and,  if  so,  can  its  structure  and/or  specific 
contents  be  identified? 

The  set  D is  frequently  taken  to  consist  of  a finite-dimensional  sub- 
space of  C(S)  generated  by  a basis  of  linearly  independent  functions 

N 

g , ...,gy,  in  which  case  D = { £ a^  g^:  a^  € Re  for  all  k}.  For  example, 

k=i 

D might  consist  of  all  polynomials  in  the  n variables  of  highest  degree  r, 
or,  in  our  utillt)  context,  D might  be  based  on  con  ..itional  univariate 
utility  functions.  When  D is  a finite-dimensional  linear  subspa.  e of  C(S), 
is  convex  and  nonempty  (e.g.,  Buck,  1959,  Theorem  2). 

Buck  (1959)  illustrates  the  above  ideas  with  a specific  example  using 
n = 2 and  N ~ 5.  Let  f(x,y)  = xy  and  let  g through  g in  C(S)  be 

1 5 


* 
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respectively  the  identity  function  (g  =1),  x,  y,  x2  and  y2.  Buck  shoves 
that  d(D,f)  * 1/4  and  that  D,  ■ {af  + (1  - a)f  : 0 < a < 1}  where  . (x,y)  ■ 

t i 2 — — 1 

(x2  + y2)/2  - 1/4  and  f (x,y)  - x + y - (x2  + y2)/2  - 1/4. 

2 

When  D is  an  infinite-dimensional  but  'small'  subspace  of  C(S),  '.t 
may  be  true  that  d(D,f)  » 0 for  every  f € C(S)  although  Df  must  then  be  empty 
for  'most'  f 6 C(S).  A useful  theorem  developed  by  Bohman  and  Korovkin 
(Korovkin,  1959;  Lorentz,  1966,  p.  7)  shows  that  d(D,f)  * 0 for  all  f € C(S) 
can  sometimes  be  established  by  showing  only  that  d(D,f)  » 0 for  a finite 
number  of  specific  functions  f.  Consider,  for  example,  the  famous  Weierstrans 
approximation  theorem  (Lorentz,  1*5&,  p.  10),  which  says  that  d(D,f)  ■ 0 for 
all  f € C(S)  when  S ■ { (t  , . . . ,t  ) : 0 < t,  < 1}  and  D is  the  space  of  all 
ordinary  algebraic  polynomials  in  t^ t^.  Using  the  Bernstein  polynomials 


Bm(f;ti,...,tn) 


m 

Z 


m 

Z 

c ■ 

n 


f(~^ m)bmk  (ti),,,bmk  (tn)#  (4) 

i n 


where  ( t ^ ) » (™)t^(l  - t^)m  \ the  Bohman-Korovkin  theorem  allows  one  to 

prove  the  Weierstrass  theorem  by  showing  that  Bm(f)  ■+■  f as  m 00  for  each 
of  the  2n  + 1 functions  1,  t^,  t£,  (i  = l,...,n). 

Although  the  above  discussion  barely  scratches  the  surface  of  approxi- 
mation theory,  it  will  suffice  for  our  present  purposes.  We  now  return  to 
our  concern  with  two-attribute  utility  functions. 


Considerations  in  Utility  Approximations 

Throughout  the  remainder  of  this  paper  we  shall  take  T = X * Y =*  [ 0 , 1 ] 2 
without  loss  of  generality  since  the  closed  and  bounded  inteivals  for  the 
two  attributes  can  be  mapped  linearly  onto  [0,1],  Two  different  types  of 


approximations  of  continuous  u on  T by  (2)  will  be  considered  according  to 
whether  the  approximation  attempts  to  provide  exact  -values  of  u at  the 
points  in  T at  which  u is  assessed  directly,  (See  the  discussion  following 
(2).)  Cases  in  which  v(x,y)  ■ u(x,y)  for  certain  designated  points  in  T or 
for  certain  conditional  utility  functions  on  X or  Y are  cases  of  approxima- 
tion by  interpolation.  Some  elementary  interpolation  methods  are  presented 
in  the  ensuing  sections. 

The  proof  method  mentioned  above  for  the  Weierstrass  theorem  is  not 
based  on  interpolation  since  the  Bernstein  polynomials,  which  in  the  present 
context  can  be  written  as 

v (x , y ) - l l u(J/m,j/m)(“)xi(l  - x)m  i(“)y;*(l  - y)“  K (5) 
i-o  j-o  J 

do  not  generally  yield  vm(i/m,j/m)  ■ u(i/m,j/m).  In  fact,  as  show:  by  i-abev 
(191A),  the  Weierstrass  theorem  cainot  be  proved  using  interpolating  poly- 
nomials. Since  the  Weierstrass  theorem  shows  that  ] | v — u | | 0 as  m -*■  00 , 

(5)  provides  a specific  instance  for  (2)  in  which  the  approximation  for  u 
becomes  precise  in  the  limit.  Nevertheless,  it  appears  (Lorentz,  1966, 
p.  102)  that  the  convergence  of  vm  to  u is  ratter  slow  compared  to  the 

convergence  obtained  by  polynomials  in  D'  , where  = { E 1 a x yJ 

u i=o  j=o  J 

a^  € Re}. 

Approximations  (2)  can  also  be  classified  according  to  whether  v 
undergoes  the  same  affine  transformation  as  does  u when  u on  the  right  hand 
side  of  (2)  is  replaced  by  u3'3  = au  -t-  b,  a J 0.  Given  v(x,y)  = If^(x)g^(y) 
as  in  (2),  we  shall  let  v3^  denote  the  tunction  obtained  from  Ef^xjg  (y) 
when  every  instance  of  u in  this  expression  is  teplaced  by  u3*3.  It  is 


cib 

important  to  realiza  that  v need  not  be  equal  to  av  + b.  We  shall  say 

<ib 

that  v is  affine  preserving  at  (a,b)  if  and  only  if  v (x,y)  = av(x,y)  + b 
for  3ll  (x,y)  € T,  and  that  v is  affine  preserving  if  and  only  if  it  is 
affine  preserving  at  (a,b)  for  all  a > 0 and  all  real  b.  Since 

v*b(x,y)  - E I[.'u(i/m,j/m)  + b] - x)m_i(“)y^(l  - y)m~^ 

M i j 1 

* avm(x,y)  + b, 

the  approximation  (5)  is  affine  preserving.  The  approximations  considered 
in  ensuing  sections  are  affine  preserving  when  b ■ 0,  but  several  are  not 
generally  affine  preserving  at  (a,b)  when  b i*  0. 

We  shall  also  be  interested  in  whether  v is  monotonic  in  x and/or  y 
when  u is  monotonic  in  x and/or  y.  If  v is  monotonic  increasing  in  x whenever 
u is  monotonic  increasing  in  x,  and  if  v is  monctonic  decreasing  in  x whenever 
u is  monotonic  decreasing  in  x,  then  we  shall  say  that  v is  monotoricity 
preserving  in  x.  A similar  convention  holds  for  y.  In  addition,  v is  said 
to  be  monotonicity  preserving  if  it  is  raonotonicity  preserving  in  each 
variable. 


3.  ELEMENTARY  APPROXIMATIONS 

This  section  examines  simple  additive,  multiplicative,  and  additive- 
multiplicative  forms  for  v.  Conditions  for  > on  $ under  which  these  forms 
are  exact  (i.e.  j | v - u||  =0)  when  the  inactions  involved  in  (2)  are  properly 
aligned  are  presented  in  Fishburn  (1974)  and  will  not  be  repeated  here  with 
the  exception  of  a comment  folio  'ng  the  proof  of  Theorem  5. 
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The.  Simple  Additive  Form 

We  begin  with  the  simple  additive  form.  Recall  that  T = X x Y = [0 , 1 3 2 . 
THEOREM  1.  Given  fixed  (x  ,y  ) t [0,1] 2 suppose  that 

l l 

v(x,y)  = u(x,y  ) + u(x  ,y)  - u(x  ,y  ),  for  all  (x,y)  € T.  (6) 

i l i 

Then  v(x,y)  ■ u(x,y)  If  x = x or  y = y_,  and  v is  affine  preserving  and 

monotonicity  preserving.  In  addition,  with  W = max  u(x,y)  - min  u(x,y), 

T T 

(x  ,y  ) can  be  chosen  so  that 

i l 

(a)  | | v - u | | <_  (j/2)W, 

(b)  | |v  - u|  | < W .if  u is  monotonic  in  either  x or  y, 

(c)  | | v — u|  | £ (1/2 )W  if_  u Is  monotonlc  in  both  variables. 

If  u is  constant  then  W » 0 and  (a)  through  (c)  hold  with  equality.  In 

general  | |v  - u| | “0  iff  u(x,y)  + u(x  ,y  ) * u(x,y  ) + u(x  ,v)  for  all 

i l t 

(x,y)  € T,  regardless  of  how  (x  ,y  ) is  chosen.  It  may  also  be  noted  that 

i i 

th/  bound  ip  (c)  is  satisfied  when  the  constant  approximation  v(x,y)  = 

(max  u + min  u)/2  is  usad ,,  regardless  of  whether  u is  monotonl-:.  However, 
this  approximation  does  not  have  the  font  of  (6). 

Proof . The  assertions  in  the  first  part  of  Theorem  1 are  obvious  from 
(6).  The  heavy  line  segments  in  Figure  1(a)  show  where  v must  equal  u.:  To 

Figure  1 about  here 

prove  the  latter  part  of  the  theorem,  observe  that 

I v(x,y)  - u(x,y)|  = ju(x,y)  + u(x  ,y  ) - u(x,y  ) - u(x  ,y)|, 

: i i 

! j v - a | | = max  |v(x,y)  - u(x,y)| 

T 


* 


1C 

For  notational  convenience  suppose  that  u is  not  constant  and  let  max  u(x,y) 

T 

= 1 and  min  u(x,y)  = 0.  Choose  (x  ,y  ) so  that  u(x  ,y  ) = 1/2.  Then 

rjt  X i j.  1 

|| v - u | | 3/2.  If  u is  monotonic  in  x then,  when  (x,y)  is  in  region  I of 

Figure  1(a), 

|v(x,y)  - u(x,y) | £ max  { |u(x,y)  - u(x^ ,y) | , |u(x,y  ) - u(x,y^ ) | } <_  1, 

and  similar  calculations  in  the  other  three  regions  show  that  | |v  - u| | <_  1. 

The  same  conclusion  holds  if  u is  monotonic  in  y.  Finally,  suppose  that  u 
is  monotonic.  We  consider  explicitly  only  the  case  in  which  u increases 
in  x and  decreases  in  y:  other  cases  are  left  to  the  reader.  By  examining 

signs  of  u differences  it  follows  that,  when  (x,y)  is  in  region  I of  Figure 
1(a), 

|v(x,y)  - u(x,y)|  min  {max  (u(x  ,y)  - u(x,y),u(x  ,y  ) - u(x,y  )}, 

max  {u(x,y)  - u(x,y  ),u(x  ,y)  - u(x  ,y  )}} 

i j.  1 

s max  {1/2, min  {u(x  ,y)  - u(x,y) ,u(x,y)  - u(x,y  )}}  - 1/2 

- i i 

since  min  (u(x  ,y)  • u(x,y) ,u(x,y)  - u(x,y  )}  < 1/2.  When  (x,y)  is  in  region 

! 1 ~ 

II  of  Figure  1(a),  |v(x,y)  - u(x,y)|  <_  min  {max  {u(x,y)  - u(x  ,y),u(x,y  ) - 
u(x  ,y  )},  max  {u(x,y)  - u(x,y  ),u(x  ,y)  - u(x  ,y  )})  c max  11/2, min  {u(x,y)  - 

1 V ! 1 ‘ “ 

u(x  ,y),u(x,y)  - u(x,y  )}}  = 1/2  since  1/2  < min  {u(x  ,y) ,u(x,y  ) } ' max  {u(x  ,y), 

: i — — i 

u(x,y  )}  u(x,y)  1.  The  region  III  and  IV  analyses  are  similar  respectively 

Wb 

to  th  se  for  regions  I and  II.  The  affine  transformation  u on  u for  which 
max  u = 1 and  mm  u = 0 then  gives  conclusions  (a),  (b)  and  (c)  as  stated 
in  the  theorem.  Q.E.D. 

The  preceding  proof  chose  (x  ,y  ) so  that  u(x  ,y  ) is  midway  between 

1 

min  u and  max  u,  Despite  the  fact  that  many  points  in  T have  a u value  that 
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is  midway  between  the  extremes  of  u,  it  does  not  follow  that  one  of  these 

will  minimize  | |v  - u| | when  v is  given  by  (6).  Although  it  is  not  generally 

possible  to  select  (x  ,y  ) to  minimize  I I v - u I I without  knowing  u on  all  of 

li 

T,  there  are  cases  in  which  this  might  be  done  when  u is  not  completely 
known.  To  illustrate  this,  we  shall  say  that  u is  conservative  when  it  is 
strictly  increasing  in  both  variables  and 

{0  £ x < x*  £ 1,0  <_  y < y*  <_  1}  =»  u(x,y*)  + u(x*,y)  ">  u(x,y)  + u(x*,y*).  (7) 

If  approximation  v is  conservative  whenever  u is  conservative,  rhen  v will  be 
said  to  be  conservatism  preserving.  Since  (6)  gives  v(x,y*)  + v(x*,y)  « 
v(x,y)  + v(x*,y*),  it  is  not  conservatism  preserving. 

A plausible  example  of  conservatism  (Fishburn,  1973;  Richards,  1973) 
arises  when  the  pairs  in  T are  two-period  income  streams.  Then  (7)  holds  if 

the  even-chance  gamble  between  (x,y*)  and  (x*,y),  which  ensures  one  of  the 

larger  amounts  x*  or  y*,  is  preferred  to  the  even-chance  gamble  between 
(x,y)  and  (x*,y*),  which  could  (pr.  = 1/2)  result  in  the  lower  amounts  x and 
y in  both  periods. 


THEOREM  2 . Suppose  that  u is  conservative  and  v is  given  by  ( 6 ) . Let 
| i v - u | | j denote  the  value  of  | |v  - u|  | when  (x  ,y  ) is  used  as  the 
fixed  point  in  (6)  and  let 


A - u(l,0)  + u(0,l)  - u(0,0)  - u(l,l). 


Then  (x  ,y  ) can  be  chosen  so  that  A/4  <_  | |v  - u 
impossible  to  have  | |v  - ul  | <"  A/4.  Also  let 


. , < A/3,  and  it  is 

(x  ,y  ) 

i i 


= (-l)i+J (u(x, j)  + u(i,y)  - u(x,y ) - u(i,j)] 


(8) 


L 


1 
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^or  d ^ {0*1}  afld  j € {0,1}.  Then  cp^  > 0 for  each  of  the  four  (i,j)  pairs , 
and 

I lv  “ ul  I (x,y)  = max  {cP00(x»y)»9u(x  y)»9oi(x,y),cpio(x,y)}.  (9) 

Moreover,  there  exists  a unique  (x*,y*)  € T that  satisfies  cp  (x,y)  * 

o o 

j <x»y>  and  1 (x,y)  = «P  (x,y) , and  the  unique  x*  and  y*  are  specified  by 

u(x*,l)  - u(x*,0)  » j [u(l,l)  - u(l,0)  ] +~  [u(0,l)  - u(0,0)  ] (10) 

u(l,y*)  - u(0,y*)  = | [u(l,l)  - u(0,l) ] + j [u(l,0)  - u(0,0)].  (11) 

Given  (x*,y*)  as  specified: 

(i)  if  cpoo(x*,y*)  = cpoi(x*,y*)  then  | |v  - al  I (x*>y*}  “ A/4; 

1L  Vp o (x*»y*)  * cP((i  (x*,y*)  then  j j v - u 1 1 is  minimized  by  a 
point  that  satisfies  cp  (x,y)  - cp  (x,y);  and  cp  (x,y)  « cp  (x,y)  * cp  (x,y) 

0 a 10  0 1 10  0 0 

and  V.  (x,y)  * cp  (x,y)  = cp  (x,y)  are  respectively  satisfied  bv  points 
0 1 10  11  — - — L 

Interior  to  T at  which  | | v - u j | < A/3; 

(iii)  i_f  cp  (x*,y*)  > cp  (x*,y*)  then  | | v — u j 1 is  minimized  bv  a 
0 0 0 1 - 

point  that  satisfies  cp  (x,y)  = cp  (x,y);  and  cp  (x,y)  = cp  (x,y)  - cp  (x,y) 

0 0 11  o C : l & 1 ' 

and  9ofl(x,y)  = cp^(x,y)  = cp^(x,y)  are  respectively  satisfied  by  points 
interior  to  T at  which  | |v  - u| | < A/3. 

The  part  of  Theorem  2 that  precedes  (8)  simply  summarizes  assertions 
spelled  out  in  greater  detail  following  (8).  Using  (7),  the  four  equations 
of  (8)  give  the  values  of  | v - u | at  the  four  corners  of  T when  (x,y)  is 
the  fixed  point  m (6).  Equation  (9)  says  that  the  largest  value  of  jv  - u| 
occurs  at  one  of  the  four  corners  of  T for  every  choice  of  fixed  point  in 
(6).  Since  as  is  easily  checked, 
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t’’ 


\ 


i 

\ 


9 (x,y)  + 9 (x,y)  + cp  (x,y)  + cp  (x,y)  = A, 

00  11  01  10 

it  is  impossible  to  have  | |v  - u| | < A/4.  The  final  part  of  the  theorem 

shows  how  a fixed  point  for  (6)  can  be  identified  so  that  | j v — u | | < A/3 

regardless  of  the  nature  of  u so  long  as  it  is  conservative.  The  interior 

point  qualification  is  used  in  (ii)  and  (iii)  since,  for  example,  cp  (0,0)  = 

o 1 

Q (0,0)  = 9^ (0,0)  = 0,  in  which  case  ||v  — u|j^o  0j“9  (0,0)  « A. 

Proof . Using  (7)  it  is  easily  seen  that  9^  _>  0 for  all  i,j  6 {0,1} 

and  that  9i^(x,y)  > 0 whenever  (x,y)  is  on  the  interior  of  T.  To  verify  (9), 

suppose  first  that  (x  ,y  ) is  the  fixed  point  for  (6)  and  that  (x  ,y  ) 

11  2 2 

lies  in  region  III  of  Figure  1(a),  with  x < x < 1 and  y < y <1,  and 

1 2 “ 12  — 

(x  ,y  ) f*  (1,1).  Then,  by  conservatism, 

2 2 

u(x  ,1)  + u(l,y  ) > u (x  ,y  ) + u(l,l) 

2 2 — 2 2 

u (x  ,1)  + u (x  ,y  ) > u (x  ,y  ) + u(x  ,1) 

1 22—12  2 

u(x  ,y  ) + u ( 1 , y ) > u (x  ,y  ) + u(l,y  ) 

2 2 1—21  2 

with  at  least  one  strict  inequality.  Addition  of  these  three  inequalities 

then  gives  u(x  ,1)  + u(l,y  ) - u(x  ,y  ) - u(l,l)  •>  u(x  ,y  ) + u(x  ,y  ) - 

iii  2 - . 2 

u(x  ,y  ) - u(x  ,y  ),  or  9 (x  ,y  ) ">  v(x  ,y  ) - u(x  ,y  ) Consequently 

1 i 2 2 1 l 2 2 2 2 

I v(x,y)  - u(x,y) | is  maximized  in  region  III  at  (x,y)  = (1,1).  Similar 

analyses  in  each  of  the  other  three  regions  of  Figure  1(a)  shows  that, 

regardless  ot  (he  choice  of  (x  ,y  ) for  (6),  |v  - u|  is  maximized  at  one  of 

1 1 

the  four  corners  of  T,  and  (9)  then  follows  immediately. 

To  verify  (10)  and  (11)  and  prepare  for  the  final  assertions  in  the 

theorem,  we  write  9 (x,y)  = 9 (x,y)  and  9 (x,y)  = 9 (x,y)  in  terms  of 

O'.  ’ 0 1 ■ 0 

u,  using  (8),  to  obtain  respectively 


14 


V 
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V 

4 

l1 

tf 


\ 


u(l,y)  - u(0,y)  - [u(l,l)  " u(0,0)]  - [u(x,l)  - u(x,0)],  (12) 

u(l,y)  - u(0,y)  = [u(l,0)  - u(0,l)]  + [u(x,l)  - u(x,0)].  (13) 

Conservatism  implies  that  u(l,y)  - u(0,y)  decreases  in  y and  that  u(l,x)  - 
u(0,x)  decreases  in  x.  Thus,  for  each  x,  (12)  will  be  satisfied  by  a unique 
y that  decreases  as  x increases,  with  y = 1 when  x = 0 and  y = 0 when  x = 1. 
Similarly,  (13)  has  y = 0 when  x = 0 and  y ■ 1 when  x » 1 with  the  unique  y 
solution  for  each  x increasing  as  x increases.  The  resulting  curves  in  T 
described  by  (12)  and  (13)  are  shown  in  Figure  2.  Their  unique  point  of 

Figure  2 about  here 

intersection  is  (x*,y*)  as  specified  in  (10)  and  (11).  This  point  is  given 

by  the  joint  s lution  of  (12)  and  (13).  If  cp  (x*,y*)  = cp  (x*,y*)  then  all 

0 0 0 1 

four  cp^  are  equal  to  A/4  at  (x*,y*)  and  this  point  uniquely  minimizes 
| | v — u | | , as  specified  in  (i)  of  the  theorem,  following  (11). 

The  following  lemma,  which  is  a variation  on  the  theme  of  the  first 
part  of  this  proof,  will  be  used  in  dealing  with  (ii)  and  (iii). 

LEMMA  1.  If  (i,1)  f (x  ,y  ) f (x  ,y  ) and  (x  ,y  ) is  in  the  rectangle 

j l 2 2 2 2 

(or  straight  line  segment  if  i ■ x or  j - y ) two  corners  of  which  are 

A u 

(x  ,y  ) and  (i,j),  then  cp , (x  ,y  ) > cp  (x  ,y  ). 
i ij  i : ij  2 2 

Proof.  Suppose  first  that  (i,i)  = (1,1)  with  (x  ,y  ) and  (x  ,y  ) as 

: 22 

pictured  in  Figure  1(a).  Then,  by  conservatism, 

u(x  ,1)  + u(x  ,y  ) > u (x  ,y  ) + u(x  ,1) 

1 2 2 — . 2 * 

u(x  ,y  ) + u(x  ,y  ) > u(x  ,y  ) + u(x  ,y  ) 

1 2 2 — : ! 2 2 

u(l,y  ) + u (x  , y ) > u (x  ,y  ) + u(l,y  ) 

2 2~2l  2 
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with  at  least  one  strict  inequality.  Addition  of  these  inequalities  gives 

u(x  ,1)  + u(l,y  ) - u(x  ,y  ) - u(l,l)  > u(x  ,1)  + u(l,y  ) - u(x  ,y  ) - 
1 111  2 2 2 2 

u(l,l),  or  cp  (x  ,y  ) > cp  (x  ,v  ).  The  other  three  regions  in  Figure  1(a) 
11111122 

are  handled  in  a similar  fashion.  Q.E.D. 


We  now  return  to  (ii)  and  ('it)  of  the  theorem.  Consider  the  point 
labeled  Q in  Figure  2.  By  Lemma  1,  cp  (Q  ) > cp  {?'  ) and  cp  (Q  ) > cp  (Q  ) , 

0 000  O'i'l  .11  1 i 0 

and,  since  cp  (Q  ) * cp  (Q  ) , cp  (Q  ) > cp  (Q  ) . Similarly,  using  Q end 
001  111000  1*0  0 

Q , cp  fQ)>cp  (Q).  By  Lemma  1,  cp  (Q  ) and  cp  (Q  ) will  be  reduced  when 
2100  010  0 ■ 0 100 

the  fixed  point  for  (6)  io  moved  from  Q in  the  direction  of  the  arrow 

emanating  from  Q , and  therefore  | |v  - u , | cannot  be  minimized  by  taking  Q 
o o 

as  the  fixed  point  for  (6).  Similar  recul , s apply  in  the  other  three 


regions  of  Figure  2.  Consequently,  min  |v  - u 


(x,y) 


must  occur  at  a point 


on  one  of  the  cp  = cp  and  cp  - cp  curves.  We  have  already  noted  that 

00  1 v 01  i 0 

the  (unique)  minimizing  point  will  be  (x*,/*)  if  cp  (x*,y*)  = cp  (x*,y*) 

0 0 0 . 

Suppose  as  in  (ii)  that  cp  (x*.y*)  < cp  Then,  by  Lemma  1, 

00  0 i 


V “ u , * .*  < 

I I I (x*,y*) 


Iv  - uMrx  y)  for  every  (x,.  . * (x*,y*)  on  the  cp^  =*  cp^ 


curve.  (A  move  from  (x*,y*)  towards  (0,1)  will  inctease  cp  , which  was  one 


of  the  maximizing  cp^  at  (x*,y*).)  Therefore  a fixed  ptint  for  (6)  that 
minimize.'  I |v  - u I I must  lie  on  the  cp  = cp  curve,  by  L ’pothesis  in  this 

OilO 

paragraph,  cp  (x*,y*)  = cp  (x*,v*)  < cp  (x*,y*)  ® cp  (x*,y*).  As  the  fixed 


oo  1 1 oi  i ; 

point  for  (6)  moves  upwards  from  (x*,y*)  to  (1,1)  along  the  cp  = cp  curve, 

0 1 0 

cp  (x,y)  increases  continuously  up  to  A = cp  (1,1),  cp  (x,y)  dec r vases  to 
0 0 0 0 1 1 

0 * cp  (1,1),  and  cp  (x,y)  and  cp  (x,y),  which  are  equal,  may  fiuct’  ate  but 

. i o : 10 

eventually  arrive  at  zero  at  (1,1),  It  follows  from  continuity  that  ’here 

is  a point  on  this  part  of  the  cp  = cp  curve  where  cp  = cp  = cp  > ' 1 4 > 

0 1 c o '0  oo 

cp  >0,  e.g.  Q , an!  with  this  point  as  the  fixed  pcint  for  (6)  we  obta-> 
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| ! v — u | | < A/3  with  the  use  of  (9).  The  remaining  parts  of  the  proofs  of 
(ii)  and  (iii)  are  similar,  Q.E.D. 

If  (7)  is  changed  by  revers1'  % the  inequality  to  u(x,y)  + u(x*,y*)  > 
u(x,y*)  4-  u(x*,y),  which  might  suggest  a 'daring'  u instead  of  a 'conservative' 
u,  then  an  obvious  correspondent  to  Theorem  2 follows  under  appropriate  sign 
changes . 

Theorem  2 illustrates  typical  concerns  of  approximation  theory  as 

outlined  after  (3).  For  Theorem  2,  D is  the  subset  of  C(T)  whose  functions 

are  given  by  (6)  as  the  fixed  point  (x  ,y  ) ranges  over  T.  Unlike  most 

typical  cases,  D depends  explicitly  on  u.  Theorem  2 shows  that  A/4  _<  d(D,u) 

< A/3,  and  (9)  implies  that  is  nonempty.  Functions  in  or  near  to  those 

in  D were  identified  in  the  latter  part  of  the  theorem, 
u 

The  S-mple  Multiplicative  Form 

The  basic  multiplicative  approximation  for  u on  T - [0,1] 2 can  be 
expressed  as  v(x,y)  = f(x)g(y),  as  in  (2).  Although  it  is  not  necessary 
to  align  v with  u In  any  specific  way,  we  shall  consider  the  case  in  which 
v(x,y)  = u(x,y)  whenever  x = x or  y = y , where  (x  ,y  ) Is  a fixed  point 

! ' i 

in  T.  This  coincides  with  our  alignment  of  the  additive  approximation  (*> ) < 
Given  v(x,y)  * f(x)g(y)  with  u scaled  so  that  u(x^,y  ) t 0,  the  specified 
alignment  implies  that 

u(x,y  )u(x  , y) 

v(x,y)  = u(-^  ~y"j , for  all  (x,y)  € T.  (14) 

■ i 

Although  this  looks  quite  different  than  (6)  and  indeed  is  in  most  cases, 
we  s^all  see  momentarily  that  (6)  is  a limiting  case  of  (14). 
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If  u has  constant  sign  then  v as  given  by  (14)  is  monotonia  in  x or  y 
when  u is  monotonic  in  x or  y,  but  if  0 is  in  the  interior  of  the  image  of 
u then  monotonicity  preservation  does  not  generally  hold.  Another  difference 
between  (14)  and  (6)  is  that  (14)  is  not  generally  affine  preserving.  In 
particular,  when  (14)  is  used  and  au(x  ,y  ) + b 4 0, 

l i 


ab.  . 
v (x,y) 


[au(x,y  ) + b] [au(x  ,y)  + b] 

1 1 

au(x  ,y  ) + b 

l l 


= av(x,y)  + b - 


ab[u  - uu,y  )][u  - u(x  ,y) ] 


J_L 


u (au  + b) 

u ii 


as) 


where  for  convenience  we  define  u = u(x  ,y  ).  This  shows  that  (14)  is 

u l l 

affine  preserving  at  (a,b)  if  b = 0,  but  it  is  affine  preserving  when  b 4 0 

only  under  very  special  conditions,  i.e.  when  u = u(x,y  ) or  u = u(x  ,y) 

u i 1 1 i 

for  all  (x,y). 

The  essential  nature  of  v remains  unchanged  if  it  (rather  than  u) 
undergoes  a positive  affine  transformation.  In  particular,  if  we  subtract 
b from  both  sides  cf  (15)  and  then  divide  by  a,  we  obtain 


ab , , 

iJmI. 


v(x,y)  + { 


u (au  + b) 


->l 


1 1 


1 1 


- u(x,y  )][u 
i i 


- u(x  ,y ) ] . (16) 


For  convenience,  let  K = -b/[u  (au  + b)]  and  let  v,.(x,y)  denote  the  left 

i • i K 

hand  side  of  (16).  Then  (16)  can  be  written  as 


u(x,y  ) u (x  ,y) 

v„(x,y)  = J—.—1 + K[u  - u(x,y  )][u  - u(x  ,y)],  (17) 

r.  u , • i , , 

where  K is  any  real  number  other  than  -1/u  < The  latter  value  of  K is 

i * 

forbidden  since  it  corresponds  to  b = or  b = -00  in  (16).  Any  other  real 
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value  of  K is  obtainable  from  K = -b/[u  (au  + b)],  and  infinite  K is 

n 11 

forbidden  by  the  proscription  against  au  + b - 0 in  writing  (15).  It  is 

1 1 

easily  checked  that  vT,(x,y)  = u(x,y  ) + u(x  ,y)  - u(x  ,y  ),  which  is  the 

K li  li 

additive  form  (6),  when  K in  (17)  is  set  equal  to  the  forbidden  value  of 
-1/u  . Hence,  by  choosing  K for  (17)  arbitrarily  near  to  -1/u  , the 

tr  1 11 

multiplicative  approximation  (17)  becomes  arbitrarily  close  to  the  addit ive 

approximation  (6);  the  convergence  of  (17)  to  (6)  as  K approaches  -1/u  : 

1 

easily  seen  to  be  uniform.  Heme  all  results  stated  for  the  additive  case 
apply,  in  the  limit,  to  the  multiplicative  case. 

It  is  important  to  note  that  the  utility  function  u as  used  in  (15), 
(16)  and  (17)  is  precisely  the  same  function  used  in  (14).  Equation  (1?) 
simply  describes  the  family  of  all  basic  multiplicative  approximations — 
uniq  ui  up  to  isomorphism  under  positive  affine  transformations  on  the 
approximations — that  correspond  to  different  ways  of  choosing  an  ciigir  and 
scale  unit  for  u and  that  render  the  approximation  exact  when  x * x c:  y 
y . The  parameter  K in  (17),  unrestricted  except  by  K ^ -1/u  , describes 

l - 

the  different  approximations  in  this  family.  Naturally,  if  the  fixed  pm: 
(x  ,y  ) used  in  (14)  or  (17)  is  changed,  then  a different  family  of  mulr.- 
plicative  approximations  is  described  by  (17). 

Thu  multiplicative  approximation  is  more  flexible  than  the  additive 
approximation  in  the  sense  that,  in  addition  to  x and  y , it  has  the 

) J 

parameter  K that  can  be  manipulated  in  fitting  a simple  multiplicative 
approximation  to  u.  We  shall  investigate  aspects  of  this  flexibility  in  t he 
next  several  theorems.  The  first  of  these  follows  immediately  from  -be 
discussion  tollowing  (17).  Recall  that  u = u(x  ,y  )< 
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THEOREM  3.  Suppose  u is  naturally  additive,  so  that  | | v — u | j “0 

when  v is  specified  by  (6) , and  u^  ^ 0.  Then,  for  every  5 > 0,  there  is 

a K i -1/u  such  that  | |vt,  - u|  | < 6,  where  vt,  is  specified  by  (17). 

1 1 «■  »■  ————————— 

In  other  words,  any  additive  (and  continuous)  utility  function  u on  T 
can  be  approximated  arbitrarily  closely  by  an  appropriate  multiplicative 
function.  Note,  however,  that  the  converse  of  this  is  not  true.  Consider, 
for  example,  the  naturally  multiplicative  function  u(x,y)  - xy.  The  general 
form  of  (6)  for  this  case  is  v(x,y)  * cx  + dy  - cd,  with  c,d  € (0,1].  As  in 
Buck's  example  of  the  preceding  section,  the  smallest  value  of  | | v — u | | 
obtainable  in  this  case  is  1/A,  which  occurs  when  c “ d =»  1/2. 

Another  indication  of  the  flexibility  permitted  by  K is  given  by  the 
following  theorem,  which  discusses  the  possibility  of  making  the  multiplica- 
tive approximation  exact  at  a point  (x  ,y  ) for  which  x ^ x and  y + y , 

2 2 2 1 2 1 

as  pictured  in  Figure  1(a).  For  notational  convenience  we  extend  the  previous 

convention  of  writing  u(x  ,y  ) as  u by  defining  u. . according  to 

li  u ij 

u^  = u(xi,y^).  (18) 


THEOREM  4.  Given  u with  u ^ 0,  let  v^,  be  defined  by  (17)  and  let 
(x  ,y  ) be  a point  in  T at  which  x ^ x and  y i y . Then  v„(x  ,y  ) = 

2 2 2 J 2 J K 2 2 

u(x  ,y  ) ■ u for  some  K j -1/u  if  and  only  if  either: 

2 2 22  1 1 4 

(a)  {u  = u and  u = u } or  {u  = u and  u = u },  in 

1 1 2 1 12  22  11  12 21  22  

which  case  v„(x  ,y  ) = u fer  every  K;  or 
K.  2 2 22  

(b)  u ^ u ^ u and  u +u  ^ u +u,in  which  case 

1 2 1 1 2 1 1 1 2 2 1 2 2 1 

v„(x  ,y  ) = u if  only  if  K = [u  u - u u ]/[u  (u  - u ) 

K 2 2 2 2 1 1 2 2 1 2 2 1 1 1 1 1 1 2 

(u  - a )]. 


11  2 1 


1<  v *^r>i 
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Proof.  To  satisfy  v„(x  ,v  ) = u we  require 

K.  2 2 22 


u u + Ku  [u  -u][u  - u ] = u u , (19) 

2112  1111  21  11  12  1122 


using  (17).  If  u = u then  (19)  holds  if  and  only  if  u = u , and  if 
11  21  12  22 

u » u then  (19)  holds  if  and  only  if  u = u . In  either  case, 

1112  2122 

v..(x  ,y  ) = u for  every  K.  The  only  other  possibility  for  u is  u 
K.  2 2 22  1 1 1 2 

u and  u i u , in  which  case  (19)  holds  for  the  unique  K specified  in 
11  21  11 

alternative  (b)  of  the  theorem.  This  value  of  K equals  -1/u  if  and  only  if 

u [u  u - u u 1 = -u  (u  - u ) (u  - u ) , which  reduces  to 
111122  1221  1111  12  11  21 

u + u = u + u . Q.E.D. 

11  12  12  21 


Given  (x  ,y  ) as  the  fixed  point  in  (6)  or  (17)  with  u 0,  and  given 
ii  l 

x i x and  y i y , Theorem  4 shows  chat  neither  (6)  nor  (17)  can  be  made 
2 12  1 

exact  at  (x  ,y  ) if,  and  only  if,  either  (u  “ u and  u + u ) or 
22  2111  2212 

u =■  u and  u u ) . 

12  11  22  21 

Our  next  theorem  parallels  the  final  part  of  Theorem  1 in  describing 
'best'  upper  bounds  on  | |v  - u| | when  the  multiplicative  form  (14)  is  used. 


Instead  of  using  (17)  we  shall  work  directly  with  (14),  taking  min  u(x,y)  = r 


and  max  u(x,y)  = r + 1,  and  present  the  bounds  as  functions  of  r.  The  choice 


of  r corresponds  to  the  choice  of  K in  (17)  when,  for  example,  u is  fixed  with 
min  u = 0 and  max  u = 1.  Because  the  midpoint  r + | of  the  interval  for  u 
equals  0 when  r = -1/2,  we  consider  only  r > -1/2  explicitly. 


THEOREM  5.  Suppose  v is  specified  by  (14)  with  min  u(x,y)  = r, 
max  u(x,y)  = r + 1,  and  u ^ 0.  Then: 

i i 

3 y.2  ^ 3r  ^ 

(a)  If  r ^ C,  it  is  always  possible  to  have  | jv  - u|  j <_  -~T'  + 2r  + p 
and  the  value  of  u = u(x  ,y  ) that  assures  this  bound  is  u = (2r2  + 2r  +•  1)/ 

ii  ii i i 

(2r  + 1) ; it  -1/2  < r < 0 then  | i v - u | | < 1 can  be  assured  by  taking 

u = r + 1; 

1 1 


^ s 

He  j 

■tdj 

(b)  If  u is  monotonic  in  one  variable  and  r _>  0 then  it  is  always 
possible  to  have  |jv  - u| | < 1,  and  this  bound  is  assured  by  every  u € 

[ (r  + l)2/(r  + 2),r  + 1];  if  u is  monotonic  in  one  variable  and  -1/2  < r < 0 


< 1 can  be  assured  by  taking  u 


r + 1 


possible  to  have  [ | v — u j | <_  2r~+'  i’  anc*  this  bound  is  assured  by  taking 
u “ r + y;  _if  u is  monotonic  and  -1/2  <_  r < 0,  then  | |v  - u|  | £ 1 is  assured 
by  taking  u^  * r + 1. 

This  theorem  and  its  ensuing  proof  show  that  monotonicity  has  no  affect 
on  the  best  general  upper  bound  on  | |v  - u|  J when  the  image  of  u contains 
the  origin  and  v is  specified  by  (14).  Moreover,  the  bounds  in  Theorem  5 
are  respectively  smaller  than,  identical  to,  and  larger  than  the  additive 
approximation  bounds  from  Theorem  1 for  the  three  cases  (a),  (b)  and  (c). 

This  is  illustrated  by  Figure  3.  As  r increases,  the  bounds  for  multiplicative 


Figure  3 about  here 
v approach  the  bounds  for  additive  v. 

Proof . Throughout  this  proof  we  write  E = |v  - u|,  or 


u(x,y  ) u (x  ,y) 

E = | 5 u(x,y)  | . 


Given  r >_  -1/2,  our  objective  is  to  identify  a value  of  u between  r and 
r + 1 that  ensures  a 'best'  upper  bound  on  the  value  of  E regardless  of  the 
nature  of  u so  long  as  min  u = r,  max  u - r + 1,  and  u satisfies  the 
monotonicity  conditions  (if  any)  that  are  specified. 
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If  u is  subject  to  no  monotonicity  restrictions  and  r > 0 then 


max  E £ max 
T 


{tr  + l)2 

u 

1 1 


r,(r  + 1) 


The  right  hand  side  is  minimized  when  (r  + l)2/u  - r = (r  + 1)  - r2/u  , or 

lx  u 

when  u = (2r2  + 2r  + l)/(2r  + 1),  and  this  choice  of  u gives 
ii  u 


3r2 

max  E < -r—jr 
— 2rz 


+ 3r  + 1 
+ 2r  + 1* 


as  specified  in  Theorem  5(a).  Continuing  without  monotonicity,  suppose 

-1/2  <r<0.  If  u >0  then 
— 1 1 

„ „ f (r  + l)2  ir(r  +1)  , . , N n 

max  E < max  P L r,  — i u - (r  + 1)  }. 

— u 1 u 1 

ll  li 

The  right  hand  side  of  this  is  minimized  when  u is  maximized  at  u = r + 1, 

u H 

so  that  max  E < 1,  given  u >0.  If  u <0  then 
~ ll  n 

„ „ r i (r  + l)2  , , r(r  + 1)  , 

max  E < max  { P - (r  + 1)  i - r) ; 

— 1 u 1 u 

11  11 

and,  with  u = r,  max  E < max  {-(r  + l)/r,l}  = -(r  + l)/r.  Since  -(r  + 1)/; 

1 1 — 

1 when  -1/2  £ r < 0,  the  best  general  result  is  achieved  in  this  case  by 

setting  u = r + 1,  with  max  E * 1.  This  completes  the  prooi  of  Theorem  5(a) 
i i — 

For  Theorem  5(b)  assume  without  loss  in  generality  that  u is  increasing 
in  x.  (If  u decreases  in  x,  a change  in  variable  from  x to  1 - x leads  to 
the  same  results.)  An  analysis  of  E according  to  the  region  of  Figure  1(a) 
that  contains  (x,y)  yields  the  following  conclusions  when  r > 0:' 
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u (r  + 1) 

I.  u(x,y)  < u(x  ,y),  u(x,y  ) < u . max  E < max  ^ 

— 1 1 — 12  — U 

1 1 

j£*Cx-Ay-3-  _ u(x,y)  |}  * max  {l,(r  + 1)  (1  - r/u  )}  = 1; 
u H 

1 1 

II.  u(x  ,y)  < u(x,y),  u < u(x,y  ).  max  E < max  {l,(r  + 1) 
l — u—i  — 

l(r  + l)/u  - 1]} 

1 1 

III.  u(x  ,y)  < u(x,y),  u < u(x,y  ).  Same  as  II; 
l — u—i 

IV.  u(x,y)  < u(x  ,y),  u(x,y  ) < u . Same  as  I. 

— i l—ii 


Therefore  max  E < max  {l,(r  + l)[(r  + l)/u  -1]}*1  whenever  (r  + 1) 

— 1 1 

[(r  + l)/u  - 1]  < 1,  i.e.  whenever  u > (r  + 1) 2 / (r  + 2) , given  r > 0. 
u — li  — 

Suppose  next  that  -1/2  < r < 0.  Since  this  case  is  intermediate  between 


(a)  and  (c)  and,  in  each  of  these,  max  E < 1 is  obtained  by  taking  u 


r + 1 (see  below  for  (c)),  the  same  result  must  hold  for  case  (b) . 

For  definiteness  in  (c)  assume  that  u increases  in  both  variables. 
Given  r > 0,  an  analysis  of  E by  the  regions  of  Figure  1(a)  yields  the 


following: 


u u 

I.  u(x,y)  < {u(x  ,y),u(x,y  )}  < u . max  £ < max  {-*- *-*-  - r, 

— i l—u  — u 


iu(x,y,; 


r<u(x,y)<u 


- u(x,y) | } = u - r; 

1 1 


II.  u(x  ,y)  < {u(x,y),u  } < u(x,y  ).  max  E <•  max  {(r  + 1)(1  - r/u  ), 

i - u—i  — u 

r + 1 - u },  where  the  terms  in  braces  are  computed  using  u(x,y) 

1 1 

= u(x  ,y)  and  u(x,y)  = u(x,y  ); 
i l 

III.  u < {a(x,y  ) ,u(/  ,y)}  < u(x,y).  max  E < max  { ! u — (r  +1)1* 

; ; ~ i l — — li 


max  !--»■  ^ u(x,y)j}  - max  {r  + 1 - u , — 

u <_u(x,y)_fr+l  uu  ‘ 1 '* 

(r  -t  1 - u )}  = r (r  + l - u ;; 

* l a li 

1 1 

LV.  u(x,y  ) <■  {u(x,y),u  } < u(x  ,y).  Same  result  as  II, 

i A 1 1 


r + 1 
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Therefore,  given  r > 0, 


max  E < max  {u  - r,r  + ^(u  - r),r  + 1 - u ,—  — "(r  + 1 - u )} 

— xi  u xx  l x u ]1 

xx  xx 

■ max  + ^(u  - r) ,— — -(r  + 1 - u )} 

u 1 1 u n 

XX  IX 

* (r  + 1)  max  {1 — - l}. 

u u 
IX  IX 


The  final  expression  here  is  minimized  when  the  two  terms  in  braces  are 


equal,  i.e.  when  u ■ r + 1/2,  in  which  case  max  E < (r  + l)/(2r  +1).  To 
x x ~ 

complete  the  proof  of  (c),  assume  that  -1/2  < r < 0.  Given  u > 0 a 


regional  analysis  of  E shows  that  max  E is  governed  by  (r  + 1)  (u  ^ - r)/u  , 

whose  minimum  value  equals  1 when  u ■ r + 1.  A lower  bour<’  n E cannot 

i x 

be  obtained  by  taking  u <0.  Hence  max  E < 1 with  u ■ r + 1 when 

XI  — XX 

-1/2  < r < 0.  Q.E.D. 


A final  comment  on  the  simple  multiplicative  form  is  in  order  before 


we  consider  another  type  of  approximation.  Suppose  that  X is  utility 


independent  of  Y and  Y is  utility  independent  of  X in  the  generalized  sense 


discussed  by  Fishburn  (1974)  and  Fishburn  and  Keeney  (1974).  Then,  with 


(x  ,y  ) , (x  ,y  ) 6 T such  that  u 4 u and  u 4 u , where  u * u(x  .,y  ) 

oo, : i xi  01  liio  oi  o: 

and  u = u(x  ,y  ) as  in  (18),  it  can  be  shown  that  | jv  - u| | * 0 when  v 
TO  10 

is  defined  by 


(u  u - u u ) [u(x,y  ) + u(x  ,y)  - u ] 
iooi  ooii  i l ii 

+ (u  + u -u  -u  )u(x,y  )u(x  ,y) 
vU,y) “ U £J LS i i . (20) 

(u  - u ) (u  - u ) 

i . 0 1 » 1 10 


If  u is  naturally  additive  then  u + u = u + u and  (20)  reduces  to 

00  11  O'  10 

(6).  If  v + u 4-  u + u then  (20)  is  identical  to  (17)  when 

00  11  01  i 0 
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K=(u  u -u  u )/[u  (u  -u  )(u  -u  )],  and  this  value  of  K is 
11  oo  oi  io  ii  n 10  li  oi 

not  equal  to  -1/u  . This  shows  that,  when  u is  not  additive,  the  multi- 

1 1 

plicative  form  (14)  or  its  counterpart  (17)  corresponds  to  the  generalized 


version  of  Keeney's  notion  of  utility  independence  in  each  direction. 


The  Simple  Additive-Multiplicative  Form 


The  simple  additive  and  multiplicative  approximations  are  exact  along 


the  two  heavy  line  segments  shown  in  Figure  1(a)  when  (x^,y  ) is  the  fixed 
point  used  in  the  approximations.  We  shall  now  examine  a mixed  additive- 


multiplicative  approximation  that  is  based  on  two  fixed  points,  (x^,y^)  and 

(x  ,y  ),  and  that  is  exact  along  the  four  heavy  line  segments  of  Figure  1(b). 
2 2 

A different  approximation  that  is  also  exact  when  x 6 {x  ,x  } or  y £ (y  ,y  } 

12  12 

will  be  mentioned  rater  in  section  5. 


THEOREM  6 . Suppose  that 


v(x,y)  - f (x)  + g (y)  + f (x)g  (y),  for  all  (x,y)  6 [0 , 1 ] 2 , (21) 

1 2 3 3 


that  v(x,y)  **  u(x,y)  when  x ( (x  ,x  } or  y ( (y  ,y  ),  and  u + u ^ u + 

12  12  ll  22  1 2 

u . Then 
2 1 


v(x,y)  = — — — ^ — — — {u (x  ,y)u(x,y  ) + u(x  .yMx.y^)  (22) 

11221221 

- u(x  ,y)u(x,y  ) - u(x  ,y)u(x,y  ) 

12  2 1 


+ 

u(x  ,y)[u 

u ] 

+ u (x  ,y) [u 

- U ] 

1 22 

_ l 

2 1 i 

12 

+ 

u(x,y  )[u 

U ] 

+ u(x,y  ) [u 

• u ] 

1 22 

! 2 

2 ) 1 

2 l 

+ 

u u - U 

U } 

12  21  11 

22 

or,  equivalently, 


v(x,y) 


u(x,y  ) + u(x  ,y)  - u 
i i n 


[u(x,y  ) - u(x,y  ) + u - u ][u(x  ,y)  - u(x  ,y)  + u 


u + u - u - u 
11  22  12  21 


and  v is  affine  preserving. 


Proof.  Given  (21)  and  the  other  initial  conditions  of  the  theorem, 

substitution  of  x for  x and  then  x for  x in  (21)  gives 
1 2 

g (y)  + f (x  )g  (y)  - u(x  ,y)  - f (x  ), 

2 3 13  1 11 

8,  (y)  + f (X  )g  * u(x  ,y)  - f (x  ). 

2 3 2 3 2 1 2 

Simultaneous  solution  of  these  equations  for  the  'unknowns'  g (y)  and  g (y) 

2 3 

gives 


f (x  )[u(x  ,>)  - f (x  )]  - f (x  ) [u(x  ,y)  - f (x  )] 


83(y) 


[u(x  ,y)  - f (x  )]  - [u(x  ,y)  - f (x  ) ] 


f (x  ) - f (x  ) 
3 2 3 1 


Substitution  of  y for  y and  then  y for  y in  (21)  leads,  in  similar  fashion, 
1 2 ’ 


f (x) 

1 


f (x) 

3 


g (y  )[u(x,y  ) - g (y  )]  - g (y  )[u(x,y  ) 

— i~_2 _1 .2.  1 3 1 2 

g (y  ) - g (y  ) 

3 2 3 1 

(u(x,y  ) - g (y  )]  - [u(x,y  ) - g (y  ) ] 

1 2 Z_ L 2.  1 

g (y  ) - g (y  ) 

3 2 3 1 


g (y  )] 


In  addition  to  these  expressions  for  g , g , f and  f , the  given  conditions 

2 3!  3 

require 


u 

= f (x  ) + 

g (y  ) 

+ 

f (x  )g  (y  ) 

1 1 

1 

1 

2 

1 

3 13  1 

u 

= f (X  ) + 

g (y  ) 

+ 

f (x  )g  (y  ) 

12 

1 

1 

2 2 

3 13  2 

U 

= f (x  ) + 

g (y  ) 

+ 

f (x  )g  (y  ) 

21 

1 2 

2 

1 

3 2 3 1 

U 

22 

= f (x  ) + 
1 2 

g (y  ) 
2 2 

+ 

f (x  )g  (x  ), 

3 2 3 2 

so  that  [f(x)-f(x)]fg(y)-g(y)]=u  + u -u  -u  , which 

3 2 31  32  31  11  22  12  21 

is  nonzero  by  presupposition.  Hence  the  denominators  of  g , g , f and  f 

2 3 1 3 

do  not  vanish.  By  substituting  the  foregoing  solutions  for  g , g , f and 

2 3 1 

f into  the  right  hand  side  of  (21)  and  then  using  the  displayed  equations 

3 

for  u , u , u and  u , it  is  readily  verified  that  v is  given  by  (22). 
11  12  21  22 

The  equivalence  between  (22)  and  (23)  is  most  easily  established  by  showing 
that  the  right  hand  side  of  (23)  ’reduces*  to  the  right  hand  side  of  (22). 
The  form  of  v given  by  (23)  shows  immediately  that  v is  affine  preserving. 
Q.E.D. 

Although  (22)  guarantees  |v  - u|  =0  for  more  line  segments  in  T than 
does  either  (6)  or  (14),  analysis  of  | | v — u | j appears  to  be  considerably 
more  difficult  for  (22)  than  for  the  simple  additive  or  multiplicative 
approximations.  Because  of  this  I shall  discuss  only  one  specific  context 
for  (22)  that  has  interesting  and  easily  derived  properties.  The  form  of 
(22)  used  in  the  ensuing  theorem  might  be  thought  of  as  a boundary  model 
or  perimeter  model  since  it  guarantees  that  [v(x,y)  - u(x,y)|  = 0 whenever 
(x,y)  lies  on  the  boundary  of  r. 

THEOREM  7.  Suppose  u is  conservative  [see  (7)]  and  v is  specified  by 
(22)  with  (x  ,y  ) = (0,0)  and  (x  ,y  ) = (1,1).  Then  v is  conservative  also 

i 1 2 2 

< A/2,  where  A = u(l,0)  + u(0,l)  - u(0,0)  - u(l,l). 


and 


v 


u 
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This  theorem  is  directly  comparable  to  Theorem  2 which  uses  the  additive 
approximation.  Theorem  7 does  not  guarantee  a lower  | |v  - u| | than  does 
Theorem  2 (where  | |v  - u| | < A/3).  However,  | |v  - u| | for  (22)  can  be 
considerably  smaller  than  A/A,  and  can  equal  zero,  whereas  | |v  - u| | must 
be  at  least  as  great  as  A/A  when  (6)  is  used  for  v.  Moreover,  the  boundary 
model  form  of  (22)  has  the  attractive  property  of  conservatism  preservation. 

Proof.  The  form  of  (22)  specified  in  Theorem  7 can  be  written  as 

v(x,y)  - — (u(0,y) [u(x,0)  - u(x,l)  + u(l,l)  - u(l,0)]  + u(l,y) [u(x,l) 

- u(x,0)  + u(0,0)  - u(0,l) ] + u(x,0) [u(l, 1)  - u(0,l) ] 

+ u(x,l)[u(0,0)  - u(l,Q) ] + u(0,l)u(l,0)  - u(0,0)u(l, 1) } . 

Therefore,  when  y'  > y, 

v(x,y")  - v(x,y)  « - ^ {[u(O.y')  - u(0,y) ] [u(x,0)  - u(x,l)  + u(l,l)  - u(l,0) ] 

+ [u(l,y")  - u(l,y) ] [u(x,l)  - u(x,0)  + u(0,0) 

- u(0, 1) ] } 

> 0 

since,  under  conservatism,  A > 0,  u(0,y')  - u(0,y)  > 0,  u(l,y")  - u(l,y)  > 0, 
and  each  of  u(x,0)  - u(x,l)  + u(l,l)  - u(l,0)  and  u(x,l)  - u(x,0)  + u(0,0)  - 
u(0,l)  is  negative  unless  x 6 {0,1}  in  which  case  one  of  these  terms  is  zero. 
Therefore  v increases  in  y.  Similarly,  v increases  in  x. 

To  verify  conservatism,  or  (7),  for  v,  suppose  x < x'  and  y < y'.  Then, 
after  cancellations  and  rearrangement,  we  obtain 

v(x,y')  + v(x',y)  - v(x,y)  - v(x",y')  = 

[utO.y")  - u(0,y)  - u(l,y")  4-  u(l,y) ] [u(x,0)  - u(x,l)  - u(x',0)  + u(x",l)] 

-A 


and  this  is  positive  sf.ice  [u(0,y')  - u(0,y)  - u(l,y-')  + u(l,y)]  > 0 and 
[u(x,0)  - u(x,l)  - u(x',0)  + u(x/> , 1) ] < 0. 

Finally,  since  both  u and  v are  conservative  and  v = u on  the  border 
of  T,  if  n < x < 1 and  0 < y < 1 then 

max  {»(0,y)  + u(x,l)  - u(0,l) ,u(x,0)  + u(l„y)  - u(l,0)}  < min  (v(x,y) ,u(x,y) } 

1 max  {v(x,y) ,u(x,y) } < min  {u(0,y)  + u(x,0)  - u(0,0) ,u(x,l)  + u(l,y)  - u(l,l)}, 

so  that 

| v(x,y)  - u(x,y) j < min  {u(x,0)  - u(x,l)  + u(0,l)  - u(0,0), 

u(x,l)  - u(x,0)  + u(l,0)  - u(l,l), 

(24) 

u(0,y)  - u(l,y)  + u(l,0)  - u(0,0), 
u(l,y)  - u(Q,y)  + u(0,l)  - u(l,l)}. 

The  four  terms  in  braces  in  (24)  are  positive  under  conservatism  and  their 
sum  equals  2A.  Hence  the  smallest  of  these  four  cannot  exceed  A/2.  There- 
fore |v(x,y)  - u(x,y)|  < A/2  for  all  (x,y).  Q.E.D. 

The  following  assertions,  whose  proofs  are  left  to  the  reader,  make 

additional  connections  with  our  previous  discussion  of  conservatism.  Let 

R(x,y)  equal  the  right  hand  side  of  (24),  with  v as  in  Theorem  7.  Then 

R(x,y)  is  uniquely  maximized  at  the  point  (x*,y*/  specified  by  (10)  and  (11), 

with  R(x*,y*)  = A/2.  R(x,y)  decreases  on  each  ray  out  from  (x*,y*)  and  is 

constant  on  the  borders  of  rectangles  whose  corners  lie  on  the  9 =9 

00  ii 

and  9 * 9 lines  of  Figure  2.  Two  such  lo'i  of  constant  R are  identified 

0 1 .0 

by  dashed  lines  on  Figure  2. 
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4.  LINEAR  INTERPOLATIONS 

The  simple  forms  of  (2)  that  were  examined  in  the  preceding  section 

have  each  f^  and  as  either  the  identity  function  or  an  exprecsion  based 

on  conditional  utility  functions  for  one  of  the  variables.  In  the  present 

section  we  shall  consider  approximations  by  linear  interpolation  in  which 

the  functions  involved  in  (2)  that  are  not  based  on  utility  values  are  more 

complex  than  the  identity  function  but  nevertheless  retain  fairly  simple 

forms.  Needless  to  say,  a vast  array  of  nonlinear  interpolation  methods 

could  be  used  to  approximate  u(x,y),  but,  with  the  exception  of  a quasi- 

linear  form  that  is  mentioned  below,  we  shall  not  go  into  these. 

The  approximations  in  this  section  are  based  on  a set  {x  ,...,x  } of 

l P 

p 2 2 values  of  X and/or  on  a sec  {y ^ yq}  of  q > 2 values  of  Y,  where 

0 ■ x <x  <...<x  -1 

1 2 P (2 

0 - y <y  <...<y  *1. 

i 7.  q 

Within  the  context  of  (25),  we  define  the  . lowing  nonnegative  piecewise 
linear  functions: 


if  xi  < 


X <•  X 


i+- 


otherwise 


(i  = 1, . . . ,p  - 1) 


/I  - a1_  (x) 


c.£(x)  =/ 


a,  (x) 
0 


if  xi_1  x ^ x^  and  i > 1 

if  x.  < x <■  x.  . and  i < p 
i — — i+ 

otherwise 


(i  = 1,.  ,p) ; 
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0 


v. 


11  yj  - y t yj+I 
otherwise 


(j  * q - 1) 


B*(y) 


(I  - e.  ,<y) 


Bj(y) 


lf  yj_j  1 y 1 yj  and  j > 1 
if  yj  1 y i yJ+1  and  j < q 
otherwise 


(j  ■ 1 q). 


The  cx£  and  B*  functions  are  continuous  on  [0,1].  For  example,  with  1 < i < p, 

a£(x)  is  zero  up  to  x^_, , increases  linearly  fiom  0 to  1 between  x^_j  and  x^ 

decreases  linearly  from  1 to  0 between  x.  and  x. . , and  is  zero  after  x, , . 

i i+i  i+i 

It  is  also  useful  to  observe  that 


a*(x)  + a*+i(x)  » 1 on  [x^x^]  for  i » 1 p - i, 

8*(y)  + &*+i(y)  “ 1 on  [yj.yj+1l  for  j B 1 - 1, 

P q 

and,  more  generally,  that  E ct*(x)  n E ($*(y)  ■ 1 for  all  x and  y. 

i° i j**j  ^ 


Four  Models 

Within  the  context  of  the  foregoing  definitions,  we  shall  consider  the 
following  four  approximations  for  u: 


v (x,y) 


v (x,y) 
? 

v (x,y) 

3 


P 

E a*(x)u(x, ,y) 
i=i 


q 

E B*(y)u(x,y.) 

j = * 3 3 

v (x,y)/2  + v (x,y)/2 
2 


for  all  (x,y)  € T 

for  all  (x,y) 
for  all  (x,y,  „ i, 
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v (x,y) 
*♦ 


+ai(x)[l  - $j(y)]ui^+1  + [1.  - a^x) 36^ (y)uj+1  ^ 
+ [1  - ai(x)][l  - Bj(y)]ui+1^+1  for  all  (x,y)  6 [x^x.^3  x 


[yj»yJ+l]»  with  i » l,...,p  - 1 
j = l» • • • »q  ~ l> 


where  = uCx^y^)  as  in.  (18).  The  points  in  T at  which  these  approxima- 
tions are  exact  are  illustrated  by  Figure  4(a):  v « u along  the  vertical 


Figure  4 about  here 


lines;  v - u along  the  horizontal  lines;  and  v 

2 3 


u at  the  pq  lattice 


points  (x  ,y  ) , (x  ,y  ) (x  ,y  ).  These  assertions  follow  immediately 

ii  12  p q 

from  the  definitions  and  the  observation  that 


xi+i  - x 


x - x. 


''■<X'y)  ' CWr^)“<*i’y)  + (V~lt7>U(‘‘1+‘,y)  “ “i  - * - X‘+‘ 


(26) 


yJ±i  ' 


y - y; 


Vx’y)  " (^17)u(x’yj)  + (7^"rl")u(x’yj^)  if  yJ  - y - yj+i’ 


(27) 


j+i 


The  determinations  of  v(x,y)  for  (x,y)  € [x^.x^.]  x [y^  »yj^.i  ] are  noted  with 
reference  to  Figure  4(b).  The  value  of  v^ (x,y)  is  the  convex  linear  combina- 
tion of  u(x. ,y)  and  u(x. , ,y)  shown  by  (26);  v (x,y)  is  the  linearly  inter- 

1 X*r  * 2 

polated  value  between  u(x,y.)  and  u(x,y  ) shown  by  (27);  v (x,y)  is  the 

J J+1  a 

average  of  the  horizontal  and  vertical  interpolations  given  by  v (x,y)  and 

i 

v (x,y)  [modification  of  v to  the  form  Av  + (1  - A)v  allows  different 

2 3 12 

emphases  to  be  placed  on  the  horizortal  and  vertical  interpolations];  and 

v (x,y)  is  given  by  a weighted  average  of  the  four  u values  at  the  corners 

u 

of  [x.,x.,  ] * [y.,y..  ]•  The  quasi-linear  approximation  v can  be  thought 

1 IT-  J jTl 
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of  as  a two-stage  linear  interpolation  process  since,  with  a = a^x)  and 

8 = 8j(y)> 

v (x,y)  = av  (x. ,y)  + (1  - a)v  (x.,  ,y) 

4 4 1 4 It  1 

= oIBuy  + (1  - e>«lij+,]  + a - + (1  - 6)u1+ljj+1) 

* 6!c,uij + a • + (1 ' B)[aul,j+> + (1  - 

* 8v  (x,y  ) + (1  - g)v  (x,y  ). 

4 J 4 J+1 

Moreover,  if  u is  linear  along  each  of  the  four  border  line  segments  of  the 
rectangle  [x. ,x. . , ] * [y . »y . . , ],  then  v * v throughout  the  rectangle,  but, 

1 lTl  J J+l  D 4 

without  such  linearity,  coincidence  of  v and  v is  assured  only  at  the  four 

3 4 

corners. 

There  are  several  major  differences  between  v and  the  other  three 

4 

approximations.  First,  v , v and  v presume  that  p,  q or  p + q conditional 

12  3 

utility  functions  are  evaluated,  whereas  v uses  only  the  u values  at  the  pq 

4 

lattice  points.  Secondly,  v , v and  v are  clearly  within  the  format  of 

1 2 3 

(2),  whereas  this  is  not  at  all  clear  for  v . I leave  it  to  the  reader  to 
show  that  v can  indeed  be  viewed  as  a special  case  of  (2),  but  it  appears 

H 

that  this  can  be  done  only  when  some  of  the  ^ and  g1  functions  ha^e  dis- 
continuities. As  shown  by  the  definitions  of  v , v and  v , and  by  previous 

12  3 

discussion,  the  f^  and  g^  functions  used  therein  are  all  continuous  on  [0,1] 
Despite  the  e differences,  all  four  approximations  possess  certain 
preservation  properties  defined  earlier. 

THEOREM  8.  Each  of  v , v , v and  v is  continuous,  affine  preserving. 

i 2 3 

monotonicity  preserving,  and  conservatism  preserving . 


w<  * i 
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Proof.  Because  of  the  similarities  among  v . v . and  v , it  will 

1 2 3 

suffice  to  consider  only  v and  v . Continuity  for  v is  obvious  and  it  is 

i 4 1 

easily  seen  to  hold  for  v . Affine  preservation  is  easily  checked.  For 

4 

monotonicity  suppose  first  that  x^  _<  x < x*  £ x^+i  and  that  u is  monotonic 
in  x.  Then,  by  (26), 

v (x*,y)  - v (x,y)  - [u(x  ,y)  - u(x  ,y)],  (28) 

i i *±+i  i* 1 j- 

so  that  v (x*,y)  - v (x,y)  has  the  same  sign  as  u(x  ,y)  - u(x. ,y).  Hence 
l l i+1  i 

v is  monotonic  in  x on  the  interval  [x^,x^+J]  in  the  same  3ense  that  u is 

monotonic  in  x on  this  interval.  Since  this  is  true  for  all  [x, v 

1 1+*-  1 

is  monotonicity  preserving  in  x.  Suppose  next  that  y,  < y < y*  < y , , and 

j — — J+1 

that  u is  monotonic  in  y.  Then,  by  (26), 


xi+i  " x 

v (x,y*)  - v (x,y)  - (-= — 

l l *i+i  xi 


)[u(xi,y*)  - u(xJ,y)] 


from  which  it  follows  that  v is  monotonicity  preserving  in  y.  With  regard 

i 

to  v , if  x.  < x < x*  < x..  and  y € [y . »y , . , ] then 

4 1 — — J j + 1 

v (x*,y)  - v^(x.y)  = (aj,(x)  - a^x*)  ] [8  (y)  (u1+1  ^ - u^) 

(29) 


+ U-  VrtXVi.J+t 


- u 


i 


i>i 


with  a.(x)  > a.(x*).  Hence  if  u is  monotonic  in  x then  v is  monotonic  in  x 
ii  H 

in  the  same  sense  as  u on  the  interval  from  x.  to  x.,  and,  since  this  is 

l 1+1  ’ 

true  of  each  such  interval,  v is  monotonicity  preserving  in  x.  The  proof 

4 

that  v is  mcrotonicity  preserving  in  y is  similar. 
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To  establish  conservatism  preservation  we  work  first  within  a rectangle 

[v * * * * xi,xi+i  1 x [yj,yj+J  with  xi  1 x < **  1 xi+1  and  yj  1 y < y*  1 y.j+1»  aad 

assume  that  u is  conservative  [see  (7)].  Using  (28), 


v (x,y*)  + v (x*,y)  - v.  (x,y)  - v(x*,y*) 
1 ill 


x*  - X 


i+i 


- *:>[u(xi'y*)  + u<xi+i>y> 


- ufr^y)  - u(xi+i,y*)]. 


Since  conservatism  for  u implies  that  the  right  hand  side  of  this  equation  is 
positive,  it  follows  that  v is  conservative  in  the  rectangle.  Using  (29), 

v (x,y*)  + v (x*,y)  - v (x.y)  - v (x*,y*)  - [a . (x)  - a,  (x*)  ] [B  (y)  - 3 (y*)] 

4 4 44  13-JJ 

+ V.,j  - uij 

* “i+i.j+il’ 

and,  since  each  term  on  the  right  hand  side  is  positive,  the  left  side  is 
positive  also.  Hence  v is  conservative  in  the  rectangle.  Therefoie  both 

4 

v and  v are  conservative  in  every  rectangle  of  the  form  [x. ,x.,  1 x 
i 4 i i+i 

[y  ,y  ].  It  then  follows  without  difficulty — by  breaking  any  rectangle 
[x,x*]  x [y,y*]  into  subrectangles  according  to  the  grid  of  Figure  4(a),  and 
adding  up  the  inequalities  implied  by  conservatism  on  each  of  the  subrectangles 
that  v and  v are  conservative  throughout  T.  Q.E.D. 

1 4 

Uniform  Norm  Considerations 

Given  (25)  and  (x,y)  6 tx^,x^+1]  * ^yj,yj+i^’  aksoIute  differences 

between  u(x,y)  and  the  approximations  defined  above  are  as  follows: 
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for  all  h € [0,1].  Considered  as  functions  of  h,  w is  the  modulus  of  con- 
tinuity for  u,  and  to  and  to  are  partial  moduli  of  continuity  for  u.  These 

1 2 

definitions  presume,  nothing  about  differentiability  or  monotonicity  for  u. 

It  is  easily  seen  that  the  moduli  are  nonnegative,  nonder.reasing  and  satisfy 


max  (id  (u,h)  ,to  (u,h)}  < to(u,h)  < to  (u,h)  + w (u,h), 
1 2 — 1 2 


Moreover,  because  u is  continuous  on  a closed  and  bounded  set,  it  follows 

readily  from  standard  results  (e.g.,  Bclzano-Weierstrass  theorem,  existence 

of  convergent  subsequences)  that  each  of  to(u,h),  to  (u,h)  and  to  (u,h) 

1 2 

approaches  zero  as  h 4 0. 


THEOREM  9.  | |v£p)  -u||  40as  p4»  for  each  k 6 {1,2, 3, 4).  In 

particular, 

| |v<p)  - u|  | < w^u.l/tp  - 1))  for  k - 1,2, 

| |v(p)  — u I t < [to  (u,l/(p  - lj)  + to  (u,l/(p  - 1 ) ) ] / 2 < to(u,l/(p  - 1)), 

3 — 1 2 — 

||v^  - U | | < to  (u , 1/  (p  - D), 

H 


and  if  each  of  to,  to  and  w is  n concave  function  of  h then 
1 2 


V<P)  - »l 

! i 

Vu'2f 

v'p)  - u| 

1 <- 

[<*> 

3 

1 <4 

v(P>  - u| 

4 

1 < 

“<u'l6T 

Proof . 

The 

initial 

r)  for  k - 1,2, 


~~TT)  (u,T7~  TT)l/2l  <u.  1 


2 (p  - 1) 


), 


r)  • 


expressions  for  |v^(x,y)  - u(x,y)|  written  earlier  and  from  the  foregoing 
comments  on  the  moduli  of  continuity.  For  example,  omitting  u from  to'  for 


notational  convenience, 


|v^''x,y)  - u(x,y)|  < a^x)  | u^.y)  - u(x,y)  | + (1  - a^x))  | u(x±+i  ,y)  - 

£ a^xjw^l  - aj(x)](x^+l  - xi>)  + (1  - a^x)) 

w^a^xXx.^  - x^) ) 

1 - a . (x)  a,  (x) 

■ a.  (x)w  ( — : — ) + (1  - a. (x))w  (- r) , 

X ip-l  1 iP-l 

so  that 

I |v(p)  - u|  I < max  [aw  (*■—-”)  + (1  - a)w  ] 

1 0<a<i  iP"1  i P - 1 

< max  [otw  (— 1 .-)  + (1  - a)w  (■— -1---)  ] 
a IP1  IP 

- wj (l/(p  - 1)). 

If  w is  concave  in  h then  aw  ((1  - a)/(p  - 1))  + (1  - a)w  (a/ (p  - 1))  < 

l i i ~ 

w ( [a(l  - a)  + (1  - a)a]/(p  - 1))  < w ( [1/2 ] / (p  - 1)).  To  verify  the 
i “ l 

concavity  result  for  v , we  note  first  tha* 

4 

|v^p)(x,y)  - u(x,y)|  £ a^xJB^y)  | u±j  - u(x,y)  | + a^xXl  - Bj  (y)  > I u±> 

- u(x,y) | + (1  - ai(x))gJ (y) |ui+1  - u(x,y)| 

+ (1  - ai(x))(l  - 6j(y))|ui+i  ^+i  - u(x,y)  j 
£ a^xJB^  (y)w(max  {l  - a^x)  ,1  - Bj  (y) ) / (p  - D) 

+ ai(x)(l  - Sj  (y))w(max  {l  - a^x) , 8j  (y) )/ (p  - D) 
+ (1  - ai(x))Bj  (y)w(max  {ai(x),l  - Bj  (y ) } / <p  - D) 
+ (1  - ai(x))(l  - B^yJMmax  (a.  (x) ,$..(’•>}/ (p  - 1 


Therefore 


| |vvp  - u|  j <_  max  [a3ui(max  {1  - ct,l  - 3) / (p  - 1)) 

4 o<ct<i 

+ a(l  - 6)<a(max  {1  - a,g}/(p  - 1)) 

+ (1  - a)gw(max  {a,l  - B)/(p  - 1)) 

+ (1  - a)(l  - 8)ca(max  {a,B}/(p  - 1))], 


Suppose  that  a <_  8 £ 1 - a.  Then,  assuming  concavity  for  w, 


,<P) 


- u 


< max  [aBwC^r)  + a(l  - 0)u>(^-£) 


a, 8 


VP  - V 

+ (1  - o)5u(if-|)  + (1  - 0(1  - B)u(— f-y)] 


< max  + (1  - cOoj(^^ — J 

~ a,S  P-1  P'1 

< MX  M((l-a)[c.  + 2B(l-  B))) 

'a, 8 P"1 

“ u(- r) 


where  the  maximizing  values  are  a * 1/4  and  8 * 1/2.  Because  of  the 
symmetry  in  a and  1 - a,  and  in  8 and  1-8,  this  suffices  to  establish  the 
final  conclusion  of  the  theorem.  Q.E.D. 


Although  Theorem  9 might  be  used  in  practical  situations  to  estimate  an 
upper  bound  on  the  maximum  difference  between  v(y,y)  and  u(x,y)  for  a given 
p that  can  be  assured  by  a linear  model,  it  says  very  little  about  good  ways 
to  choose  the  x^  and/or  y for  (25)  in  attempting  to  minimize  | |v  - u| | for 
fixed  p and/or  q. 

To  illustrate  the  latter  idea,  suppose  u increases  in  both  variables 


and  v is  used  with  p = q * 3.  Then  the  best  general  assertion  that  can  be 

H 

made  for  an  upper  bound  on  | |v  — u j | is 


40 


| |v  - u | | <_  min  [max  {u(x  ,1)  - u(0,y  ),u(l,y  ) - u(x  ,0), 

4 0<X2<1  2 222  (30) 

0<v  <1  uU.l)  “ u(x  ,y  ),u(x  ,y  ) - v(0,0)}]. 

y o 2 2 2 2 


For  convenience  set  u(0,0)  = 0 and  u(l,l)  ■ 1,  and  let  V denote  the  value  of 

the  right  hand  side  of  (30).  Then,  because  of  1 - u(x  ,y  ) and  u(x  ,y  ) - 0 

2 2 2 2 

in  (30),  V >_  1/2.  In  any  event,  it  can  be  shown  that  some  point  that 

satisfies  u(x  ,1)  - u(0,y  ) » u(l,y  ) - u(x  ,0),  or  u(x  ,0)  + u(x  ,1)  » 

2 2 2 2 2 2 

u(0,y  ) + u(l,y  ),  must  be  a minimaxing  point  for  V.  There  exists  a unique 
2 2 

(x*,y*)  that  satisfies  u(x*,0)  + u(x*,l)  - u(0,y*)  + u(l,y*)  and  u(x*,y*)  - 
22  2 2 2222 

1/2.  Consequently,  if  u(x*,l)  - u(0,y*)  < 1/2  then  V = 1/2,  and  if  u(x*,l)  - 

2 2 ~ 2 

u(0,y*)  > 1/2  then  V > 1/2.  In  the  latter  case,  (x*,y*)  may  or  may  not  be  a 
2 22' 

minimaxing  point  for  V,  depending  on  the  behavior  of  u(x  ,1)  - u(0,y  ) 

2 2 

relative  to  u(x  ,y  ) along  the  curve  through  T that  gives  the  (x  ,y  ) solu- 
2 2 2 2 

tions  to  u (x  ,0)  + u(x  ,1)  ■ u(0,y  ) + u(l,y  ).  However,  the  choice  of 
2 2 2 2 

(x*,y*)  as  the  interior  point  to  use  for  v when  p * q » 3 appears  to  be 

2 2 4 

reasonable. 


5.  EXACT  GRID  MODELS 

In  concluding  our  discussion  of  approximations  for  u on  T ■ [0,1] 2 we 

shall  consider  several  approximations  that  are  exact  (v  = u)  on  both  the 

horizontal  and  vertical  line  segments  of  a grid  on  T such  as  shown  in  Figure 

4(a).  To  focus  the  discussion  we  shall  say  that  an  approximation  v is  an 

exact  grid  model  if  v(xi»y)  * u(xi>y)  for  all  y € [0,1]  for  at  least  two 

distinct  x^  and  if  v(x,yj)  = u(x,y^.)  for  all  x 6 [0,1]  for  at  least  two 

distinct  y . . 

J 
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The  only  approximation  of  previous  sections  that  is  an  exact  grid  model 

is  the  simple  additive-multiplicative  model  discussed  at  the  end  of 

section  3.  Approximations  v and  v of  the  preceding  section  are  not  exact 

3 4 

grid  models  since  they  are  exact  only  at  the  lattice  points  or  intersection 
points  of  the  grid.  On  the  other  hand,  the  simple  forms  of  section  3 can  all 
be  adapted  to  serve  as  exact  grid  models  on  any  grid  formed  from  a finite 
number  of  horizontal  and  vertical  lines  through  T by  applying  these  forms  in 
a patchwork  or  cut-and-paste  fashion  to  different  sections  of  the  grid. 

Suppose,  for  example,  that  (25)  holds  with  P ^ 3 and  q >_  3.  Then,  provided 

that  + ui+1>j+1  - - ui+1^  i 0 for  each  (i,j)  6 {1 - 1}  * 

{l,...,q  - l},  the  simple  additive-multiplicative  form  can  be  applied  separately 
to  each  rectangle  [x. ,x,,l  x [y.,y.., ].  That  is,  (23)  with  x ,x  ,y  ,y 
replaced  by  xi»x1+1  »yj  *yj+t  is  taken  t0  hold  throughout  [x^x^]  x [y^yj+J 

for  i ■ 1 - 1 and  j = l,...,q  - 1.  The  resultant  approximation  is  exact 

on  the  grid  of  Figure  4(a). 

Alternatively,  with  0=>x  <x'<x  < x'  <,..<  x < x*  x =1 

1122  P-1  P-1  P 

and  0 = y <■  y'  < y < y'  <...<  y < y'  < y =1,  one  could  use  the 

l l 2 2 'q-i  7q-i  q 

additive  form 


v(x,y)  = u(x,y')  + u(x',y)  - u(x',y')  for  all  (x,y)  6 [x^x^)  » Ly^  »yj+1 ) 


except  that  tx^>xi+1)  is  replaced  by  [x^,x^+1]  when  i = p - 1,  and  [y^.y^) 

is  replaced  by  [y^,y^+i]  when  j = q - 1.  The  resultant  approximation  is 

exact  on  the  grid  determined  by  x',x',...,x'  and  y' ,y' , . . . ,y'  , but  it 

12  p-i  iw2  q-i 

has  one  serious  disadvantage  that  is  not  shared  by  the  patchwork  adaptation 
of  the  simple  additive-multiplicative  model  based  on  (23),  and  that  is  its 
propensity  for  discontinuities  along  the  x^  and  y lines.  A similar 
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disadvantage  arises  with  an  adaptation  of  the  simple  multiplicative  model . 
Hence,  because  of  the  analytical  difficulties  that  accompany  discontinuities, 
the  only  one  of  the  three  simple  forms  of  section  3 that  appears  to  adapt 
itself  reasonably  well  to  a patchwork  format  is  the  additive-multiplicative 
form.  The  adaptation  of  this  form  is  continuous,  and  Theorem  7,  involving 
conservatism,  is  easily  generalized  to  the  patchwork  format. 

A Generalized  Multiplicative  Form 

In  the  remainder  of  this  section  we  shall  consider  a generalization  of 
the  simple  multiplicative  form  that  is  not  a patchwork  adaptation.  It  is 
based  on  m fixed  points  for  X and  for  Y subject  to 

0<x  <x  < . . .<  x <1 

“ 1 2 (31) 

0 i y,  < y2  <•••< 

and  is  derived  directly  from  (2)  and  the  restriction  that  v = u along  each 
line  determined  by  the  2m  fixed  points  in  (31).  Its  basic  form  is  given  as 

m m 

v(x,y)  = T.  E c,.  u(x,y  )u(x,  ,y),  (32) 

i=i  j*i  3 3 

where  the  are  based  on  the  m-by-m  matrix  of  u^  values  (i,j  = 1,  . ,m) . 

When  m = 1,  c =■  1/u  and  (32)  reduces  to  the  simple  multiplicative  model 

i . 1 1 4 

(14).  When  m ~ 2,  c * u /A,  c = -u  /A,  c = -u  /A,  and  c = u /A, 

n 2 2 12  21  2.  x2  22 

where  A = u u - u u , the  determinant  of  the  u.,  matrix.  Thus  the 
.1221221  ij 

m = 2 version  of  (32)  provides  an  alternative  to  the  exact  2-by-2  grid  model 
of  the  simple  ac.ditive-multiplicative  form  so  long  as  the  u^  matrix  is 


nonsingular . 
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Approximation  (32)  seems  attractive  for  several  reasons.  It  has  a nice 
analytical  form,  is  exact  on  the  grid  determined  by  the  points  in  (31) , and 
is  continuous.  Moreover,  if  u is  differentiable  then  v is  differentiable. 
However,  it  is  neither  affine  preserving  nor  monotonicity  preserving,  and  it 
does  not  submit  easily  to  analyses  of  j | v — u | | . Although  very  little  is 
known  about  ||v  - u|j  at  the  present  time,  it  is  hoped  that  further  research 
will  determine  the  conditions  under  which  (32)  gives  a good  approximation 
to  u. 

Because  of  the  absence  of  interesting  results  on  | |v  - u| |,  I shall 
present  only  the  basic  derivation  of  (32)  from  (2)  and  note  the  effects  of 
positive  affine  transformations  on  this  approximation.  This  presentation 
parallels  the  discussion  for  the  simple  multiplicative  form  in  (14)  through 
(17).  In  the  present  setting  we  shall  presume  that  m _>  2 and,  given  u and 
(31),  let  U denote  the  m-by-m  matrix  [u^l.  Also  let  be  the  (m  - l)-by- 
(m  - 1)  matrix  obtained  from  U by  deleting  its  ith  row  and  jth  column  and, 
with  det  the  determinant  function  on  square  matrices,  define 

A = det  (U) 

Ajj  = (-l)i+:i  det  (U^)  i,j  = 1 


so  that  A^j  is  the  cofactor  of  u^ . 

THEOREM  10.  Let  u and  (31)  be  given  with  m > 2,  and  suppose  that 


m 


v(x,y)  = l f,.(x)g,.(y)  for  all  (x,y)  6 [0,1]' 
k=’  K K 


(33) 


with  v(x,y)  = u(x,y)  whenever  x € {x  ,...,xm}  or_  y ( {y  .....y^},  and  that 


U is  nonsingular  (A  ^ 0) . Then 
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m m 

v(x,y)  = E E (A44/A)u(x,y.t)u(x.l,y)  for  all  (x,y)  € [0,1]2.  (34) 

i-i  j-i  J 2 


Moreover,  if  a and  b are  real  numbers  with  a > 0.  and  if  aA  + bEm 


0,  then 


ab . 


au * % - mm 

— “ v(x,y)  + K(a,b)  [A  - E E A.,u(x,y)] 

a i-i  j-i  1 


m m 


[A  - E E A..  u(x,y  )] 
i-i  j-i  3 3 


where 


(35) 


K(a,b) 


-b 


m m 

A[aA  + b E E A 
i-l  j-i 


ij 


(36) 


Expression  (35)  is  the  appropriate  generalization  of  (16)  or  (17)  for 

m ■ 1.  As  seen  by  (36),  the  only  real  value  of  K that  cannot  be  obtained 

with  admissible  values  of  a and  b is  K = -l/CAE^  A^),  corresponding  to 

a = 0,  provided  that  E^  A^  ^ 0.  If  E^  A^  = 0 then  every  rea)  K is 

admissible.  Moreover,  if  b 4 0,  then  (35)  shows  that  va^(x,y)  - av(x,y)  + b 

if  and  only  if  [A  - E^  A^  u(x^,y)][A  - E^  A^  u(x,yj)]  ■ 0,  which  holds 

if  x € {x  , ...,x  } or  y € {y  ,...,y  } since  E.  A.,  u . . = E.  A,,  u , , = A and 
i o l m i ij  i]  ] ij  ij 

^i  Aih  uij  = 0 h ^ j and  E^  A^  u^  = 0 if  k j*  i,  but  which  cannot  be 
expected  to  hol<_  otherwise. 

We  conclude  with  an  outline  of  the  proof  of  Theorem  10  since  the 
complete  proof  is  rather  long.  Assume  henceforth  that  the  first  sentence 
of  the  theorem  applies.  Substitution  of  xi  and  then  y^  into  (33)  gives 
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m 

£ ft(x.)g.  (y)  = u(x  ,y)  i = l,...,m 
k-i  K 

m 

£ fk^My-P  = u(x,yi>  J “ 1 m* 

k*i  J J 

Let  F be  the  m-by-m  matrix  whose  entry  in  row  i and  column  k is  fk(x£)>  and 

let  G be  the  m-by-m  matrix  whose  entry  in  row  j and  column  k is  g^Cy^). 

F^  and  are  defined  similarly  to  by  deleting  rows  and  columns.  Then, 
if  F and  G are  nonsingular,  Cramer's  rule  gives 

m i+k 

g,(y)  - E u(x.,y)  (-1)  det  (F  ,)/det  (F)  k=  l,...,m 

i*i  ~A 

m 1+k 

f,(x)  » E u(x,y.)(-l)J  det  (G.,)/det  (G)  k*  l,...,m. 

k j.j  J Jk 

It  then  follows  from  the  product  and  transposition  rules  for  determinants 
(with  a prime  denoting  transposition)  that  det  (F)  det  (G)  ■ det  (F)  det  ( G ') 
- det  (FG')  = det  ([E^  ^k^Xi^8k^yj ^ ^ = c*et  ^ since  the  initial 

conditions  require  u^  = E^  ^k^Xi^8k^yj ^ ' anc*  hence  that  neither  F nor  G is 
singular  with 

m .mm  , , , m 

£ fk(x)gk(y)  * J £ £ u(x,y  )u(xity)(-l)  J E det  (Fik  G -k) • (37) 

k=i  i=i  j=i  3 k=i  J 

It  can  be  shown  that  det  (d^)  = £k  det  (F^k  G^k) , and  therefore  (33)  and 
(37)  yield  (34). 

Assume  henceforth  that  (34)  holds  and  that  aA  + bE . . A, , ? 0 with 

ij 

a * 0.  Let  * [au^  + h]  with  A3*3  = det  (U3'5)  and  A3^  the  cofactor  of 
au^j  + b in  U . Then  direct  substitution  in  (34)  gives 
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vab(x,y)  = E E (Aab/Aab) [au(x,y.)  +b][au(x  ,y)  +b]. 
i=i  j=i  J 2 1 


(38) 


Let  a denote  a generic  permutation  on  let  sgn(a)  be  the  number 

of  inversions  in  a (an  inversion  occurs  when  i < j and  ) , and  let 

tt (i , k; a)  be  the  product  of  the  m - 2 u^  for  those  j t {i,k}.  It  can  then 
be  shown  that 


, mm 

Aab  - a“  A + am_1  b E E A,  , 

i-i  j-i  2 


(39) 


.ab  m-i  . , m-2 

Aij  " a Aij  + a b 


£ (-l)sgn(a)  E 7i  (i,k;a) , 
{ k : k#  i } 


(40) 


and  that 


A E (-i)sgn(a)  E tt  ( i , k ; a ) 

{k:k^i} 


m m 


Aij ( 1 1 \h} 

K=i  n= 1 


m 


m 


( E A )(  E A..), 
h-i  lh  k-  *j 


(41) 


The  last  identity  allows  one  to  conclude  immediately  that  the  sum  of  its 
left  hand  side  over  i equals  zero  and  that  the  sum  of  its  left  hand  side 
over  j equals  zero.  This  fact  along  with  substitution  of  (40)  and  (41)  in 
(38)  gives 


.ab  ab,  . 
A v (x,y) 


= [am  1 Atj  + am  2 bE(-J)sgn^  Eir(i,k;o)  J [a2  u(x,y.. )u (Xj ,y) 

+ ab(u(x,yi)  + u(xi>y))  t b2] 

= am'  Ei.  A u(x,y  )u(x  ,y)  + ara  bE  A (u(x,y  ) + u(x  ,y)) 

J J J J J J 


+ am_  b2  E. . A. . + am  bE , . 

iJ  ij  ij 

(Zh  Aih}  (Ik  \j)]/A* 


u(x>y.)u(xi,y)[Aij  Ekh  Akh 
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Multiplication  by  A and  the  use  of  (34)  then  gives 


AAab  vab(x,y)  = a™*1  A2  v(x,y)  + a”  bA[Ey  Ay  u(x,y^)  + Ey  Ay  u(x^,y)] 
+ am  1 b2  AEy  Ay  + am  bAEy  Ay  v(x,y) 

- am  b (E±j  Ay  u(x,y^ )) (E^  Ay  u(xi>y))  + am  bA2  ••  am  bA‘ 
* bA[am  A + am  bEy  Ay]  + av(x,y)A[am  A + am  * bEy  Ay] 

- am  b [A2  - A(Ey  Ay  u(x.y^)  + Ey  Ay  u(xi,y)) 

+ (Ey  Ay  u(x,y^))(Ey  Ay  u(xi>y)) ] 


which,  since  AAab  * A[am  A + am  1 bEy  Ay]  by  (39),  yields 


vab(x,y)  = c»v(x,y) 


am  b [A  - E A u(x,y  )][A  - E A u(x  ,y)] 
,y)  + b i*1 J - 


A[am  A + am  bEy  Ay] 


Since  this  is  equivalent  to  (35),  the  outline  of  the  proof  of  Theorem  10  is 
completed. 
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max  u(x,y)  = r + 1. 
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